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MEASURES OF INDUSTRY CAPACITY! 


Edward Zabel 
Princeton University 





A survey of statistical estimates of industry capacity is made, and 
some possible uses for the measures are discussed. 











INTRODUCTION 

In a previous paper the author investigated in rather general terms the meanings and 
implications of a capacity.2 In the present study an attempt is made to explore particularly 
the concepts which underlie statistical estimates of industry capacity. 

For a sizable number of American industries measures of capacity are published in 
trade journals and in various federal government periodicals. In general, for industries 
whose operations are mechanized and whose activity is governed by durable production equip- 
ment, the journals tend to publish estimates of capacity; whereas for industries in an agri- 
cultural or services classification, such estimates are rarely made. No comprehensive 
periodic survey is made of industry capacity measures, so it is not possible, without examining 
hundreds of trade journals, to determine the number of industries which do publish figures. 
But some idea of magnitude may be surmised from published studies on capacity measures. 
In 1953 the Inter-Industry Analysis Branch of the Bureau of Mines compiled capacity estimates 
for 23 Emergency-Model (192 Industry Classification) industries. The 23 E-M industries are 
primarily in the minerals field and correspond to 92 SIC (Standard Industry Classification) 
industries. Around one third of these estimates were derived from published data on capacity, 
while the others were determined from adjusted output figures.3 Also in 1953 a study was 
prepared from published capacity data on 15 industries producing basic materials such as 
steel ingots, petroleum, and paper.4 A third study, prepared in the Office of Price Adminis- 
tration in 1941, was concerned with the capacities of 71 manufacturing industries. Approxi- 
mately two fifths of the estimates were taken from trade journal sources, while the remainder 





1The research for this paper was carried out under office of Naval Research Contract 
Noonr-27009. The author is indebted to Professor Oskar Morgenstern of Princeton University 
for his guidance and suggestions. 

2Edward Zabel, "On the Meaning and Use of a Capacity Concept,''Naval Research Logistics 
Quarterly, December 1955, pp. 237-249. 

Henry Goldstein, Pierre Crosson, and Sidney Sonenblum, Capacity Estimates for 23 E-M 
Industries, Inter-Industry Analysis Branch, Item No. 36, Bureau of Mines, 1953. 

Bert G. Hickman, "Cyclical Variations in Capacity and Capacity Utilization in Selected 
Industries,"' Conference on Research in Income and Wealth, October 9-11, 1953. 


*Manuscript received June 14, 1955. 
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were derived from informal estimates of trade association officials and federal government 
industry experts. 5 

An attempt at comprehensively reporting industry capacities was made as a special 
inquiry for the United States census of manufactures for 1921 and 1923, but after the latter 
census the inquiry was dropped. As late as the Census of Manufactures, 1947, only scanty 
data on industry capacity are to be found. The data on productive capacity presented in the 
census for 1923 were derived primarily from the following entry on the questionnaire submitted 
to each manufacturer: 

“What is your estimate of the percentage of your output compared with 

your possible output if you had such demand as to require full running time?" 
From the data obtained in response to this question and presumably from other data, the 
census authorities compiled a ratio of actual output to maximum possible output for all but 
one of the 333 industries covered in 1923.6 

The census authorities were not satisfied with the results of the inquiry, and it was 
this dissatisfaction which led to the abandonment of similar inquiries in subsequent years,7 
The Census Bureau, however, must be held at least partially responsible for the failure of the 
survey since, clearly, the question they asked was not completely unambiguous. As the replies 
indicated that such ambiguity was present, the rational procedure would seem to have been to 
revise the questionnaire rather than to drop the inquiry entirely. 





THE CHARACTERISTICS OF INDUSTRY CAPACITY MEASURES® 

The published capacity measures, generally, have important common characteristics, 
but the differences which prevail among them in terms of accuracy of data and the precise 
concept which is the basis for measurement make it mandatory for the user to be well aware of 
the peculiarities of each estimate. In common, most of the measures purport to show the maxi- 
mum amount of output which can be produced during some period of time in the industry with 
the existing durable production equipment. Since the maximum output is not only a function of 
the durable equipment, however, an answer to this question cannot be derived without specifying 
what conditions are assumed to prevail in terms of the complementary inputs, the number of 
shifts worked by the labor force, and the time, if any, allowed for the maintenance and repair 
of machinery. These assumptions vary widely from industry to industry and, consequently, the 
raw measures cannot be considered comparable without adjustments for different sets of 





5Robert A.Solo, Industrial Capacity in the United States, Office for Emergency Management 
and Office for Price Administration and Civilian Supply, June 1941. 

Census of Manufactures, 1923, pp. 1201-1207. 

The misgivings were indicated in the following letter from the Chief Statistician for 
Manufactures, U.S. Bureau of the Census, to the Editor of Food Industries: 

"Many manufacturers failed to answer the question at all, and the replies made by others 
indicated that they had failed to understand the import of the question. The purpose was to 
ascertain the per cent which the actual output formed of the maximum possible output which 
could have been produced with the equipment in place; but apparently many manufacturers 
interpreted the question to refer tothe maximum possible output which could have been produced 
with the equipment in place, together with such additional equipment as could have been installed 
in the space available. It became necessary, therefore, to do agreat deal of estimating in order 
to assemble the statistics; and, since it is the function of the Bureau of the Census to compile 
and publish reliable and trustworthy statistics, not estimates of doubtfulauthenticity, the inquiry 
was dropped." 

Quoted in Edwin G. Nourse and Associates, America's Capacity to Produce, Washington, 
D.C.: The Brookings Institution, 1934, p. 169. 

8The findings in this section are based on the author's survey of approximately 35 industry 
measures of capacity published in trade journals and federal government publications. 
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assumptions. For the bituminous coal industry, for example, the capacity figure is derived by 
dividing the actual annual output by the average number of days worked and multiplying this 
figure by 308 days. The capacity estimate is thus influenced by the operating conditions and 

the economic factors affecting the supply of labor and other inputs in each year, but it does not 
allow for strikes, breakdowns, or repairs. The data, moreover, are for operating mines only, 
though idle mines may easily be opened.? A more usual practice is to assume an unlimited 
supply of labor and raw materials, but with varying conditions as to the time use of the machin- 
ery, i.e., adjustments may or may not be made for repairs and maintenance, An example is 

the wheat flour industry where capacity is computed in the following manner. A daily capacity 
for machinery is obtained from engineering data on how many barrels of wheat flour can be 
processed over a twenty-four hour period without regard for labor and other input costs. These 
figures are reported monthly, and an average of the monthly figures is derived to determine 
average daily capacity. This average is then multiplied by the number of working days in the 
year (306 during 1923-1948, and 255 from 1949 forward) to obtain annual capacity. The mills 
covered in the data produce from 90 to 95 percent of the flour. 10 

A serious drawback of capacity statistics is that frequently the complete set of assump- 
tions which are the basis for the estimates are not explicitly described. Conditions as to the 
number of shifts in operation, number of hours per shift, working days per year, and shutdown 
time for repairs and maintenance are often not listed. Another drawback for some purposes is 
that the published data generally tell us nothing about the age distribution and the comparative 
efficiency of the production units included in the capacity totals. In the same vein, the data 
frequently do not tell us if output attributable to idle or obsolescent plants and machinery are 
included in the estimates. 

So far we have discussed assumptions which are necessary to derive capacity estimates 
for single-commodity, single-process industries. Where the industry produces multiple prod- 
ucts or where multiple production processes exist, new problems arise. For the multiple- 
commodity industry where no assumptions are made concerning the product mix, capacity 
estimates must be given in terms of a transformation function showing the various maximum 
combinations of output which can be produced or conceivably in terms of the maximum amount 
of some input item which can be handled by the fixed equipment. The latter would be an accept- 
able procedure for most uses only if the maximum input flow does not vary greatly with dif- 
ferent product mixes. In petroleum refining, for example, the capacity figure for the distil- 
lation process is given in terms of the number of barrels of crude oil run to stills. This 
appears to be a reasonable method, since in distillation the product mix is very narrowly 
determined by the characteristics of the crude oil.!11 The method used in the paper industry 
is to base the capacity estimate on what is called normal expectations as to the grade and 
weight of paper desired by the market.12 For multiple-process industries various methods 
are used. An example is the cement industry where production essentially involves two proc- 
esses; (a) the production of clinker by baking kilns and (b) the conversion of this clinker into 
finished cement by grinding mills. Capacities for each process can be assigned to any one 





9Minerals Yearbook, 1949, p. 282. However, in Bituminous Coal and Lignite in 1950, 
Mineral Market Report No. 2032, Department of Commerce, p. 7, the definition of capacity is 
based on operations for 280 days. This latter definition was supplied by the American Institute 
of Mining and Metallurgy. 

10Facts for Industry Series, M-16A, U.S. Bureau of the Census. 
lMinerals Yearbook, various issues. 
2American Paper and Pulp Association, Statistics of the Paper Industry, 1951, p. 53. 
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plant, and the ratio of these different capacities may vary within wide limits from plant to 
plant. The cement industry, however, publishes official capacity figures based on only its 
annual capacity to produce clinker, 13 

Industry measures are generally designated as technological measures of capacity 
because of their noncost nature, but it is customary to classify these concepts even more 
finely. Three types of technological measure are frequently referred to: rated capacity, 
theoretical capacity, and practical capacity. There are no standard definitions for these 
terms, but in general the distinctions are the following. Rated capacity refers to an engineer- 
ing rating (or an estimate based on past behavior) of the output potential of machinery during 
some standard period of time, perhaps an hour, a day, a month, or a year. It assumes an 
unlimited and uninterrupted flow of labor and raw materials, and, if it is computed for a short 
period of time, it does not usually adjust for delays due to maintenance and repair. The rated 
capacity of a plant is generally the basic datum from which a theoretical or practical capacity 
is calculated. Theoretical capacity refers to the maximum output which can be produced 
during some period of time, usually a year, with the existing equipment, assuming uninter- 
rupted operation. For example, if the rated capacity is the amount of output which can be 
produced during a twenty-four hour period, the annual theoretical capacity for the plant would 
be this figure times the number of working days in the year. While the theoretical capacity 
thus makes no adjustments for interruptions, the practical capacity adjusts the former estimate 
for such factors as the maintenance and repair time and the usual number of shifts worked in 
the industry. The differences in meaning noted between theoretical and practical capacity hold 
fairly consistently among industries, while rated capacity may coincide with either of the for- 
mer measures. For example, when rated capacity is computed on the basis of a year it usually 
coincides with theoretical capacity, but occasionally it is adjusted for necessary interruptions 
so that it is equivalent to a practical capacity. When using capacity data, one must be alert for 
such nuances but, as we have indicated, it is not.always possible to discover the basis for the 
estimates. 

For users of capacity data another caution should be noted relating to the accuracy of 
the data. While it is difficult to make generalizations, the author encourages an awareness that 
capacity estimates are not the type of data which generally have the sanction of the major sta- 
tistical reporting agencies, such as the Bureau of the Census and the Bureau of Labor Statis- 
tics (note the comment in footnote 7). The data usually appear originally in trade journals, 
and it is difficult to ascertain what methods were used to gather the data and what importance 
they assume in particular industries. Confidence in the statistics must depend on a knowledge 
of the reporting procedure and what, if any, strategical value the statistics have in each 
industry. 14 





13Minerals Yearbook, various issues. 
4In reporting on the statistics in 1934, Edwin G. Nourse made the following statement which 
is probably equally as valid today: 

"Unfortunately many industries keep no records of capacity. Even for those that do, the 
figures that pass current in the trade require the closest scrutiny. This is particularly true of 
statements of capacity intrade directories or similar commercial publications where the report- 
ing company knows that its statement will be published. Under these conditions the return may 
be influenced by the fact that it will be read by competitors, by tax assessors, and, not least 
important, by prospective customers. Sometimes the company is tempted to understate, more 
often perhaps to overstate, its actual performance. In the case of capacity records collected 
as an incident to the allocation of markets or the establishment of production quotas, the tempta- 
tion to make the largest showing possible is strong; and for this reason the capacities reported 
by the internationalcartels and by certain American trade groups must be scrutinized with care. 

(Continued) 
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In studies on capacity, a question occasionally asked is whether one can add the capaci- 
ties of the various plants in an industry to get a meaningful industry measure or whether one 
can add industry measures to get meaningful figures.15 We have already implicitly answered 
this question in describing the concepts which are the bases for capacity measures. In a word, 
the answer depends on the use to which the data are being put. If the objective is to determine 
the maximum amount of output which can be produced in an industry under specified conditions 
as to costs and price, then, clearly, adding technical measures of plant capacity, which depend 
on an unlimited supply of labor and raw materials, will not provide a meaningful answer. In 
the same sense, if the objective is to determine the capacity output for the economy, adding 
together technological measures of industry capacities will not give a meaningful answer. On 
the other hand, if the problem is to determine an indicator of the size of an industry or the 
growth of an industry over time, the objections to adding technical plant capacities lose most 
of their force. In the same sense, if one wished to compare differential rates of growth for 
different groups of industry, comparing summed groups of industry capacity measures may be 
useful as long as the nature of the capacity measures is kept clearly in mind. The uses and 
limitations of capacity measures are discussed in more detail in a later section. 

In order to make our discussion on industry measures of capacity more pointed, in the 
next section we shall describe the estimating procedures used in some important industries. 


ESTIMATING PROCEDURES USED IN SOME BASIC INDUSTRIES 
1, The Steel Industry 

Statistics on the steel industry are originally compiled, in the main, by the American 
Iron and Steel Institute and published initially in various of the Institute's periodicals, such as 
Steel Facts and the Annual Statistical Report. Data on the capacity of the steel industry are 
published annually and are obtained by the Institute from questionnaires sent every year to each 
firm in the industry. A firm is requested to make individual reports for every plant that it 
owns, and the coverage includes all plants which are primarily engaged in producing pig iron 
and ingots for casting. 

The total annual steel- making capacity in the United States is announced once a year by 
the Institute, generally as of January 1st, and sometimes at mid-year as well. This practice 
was begun at the start of 1914, on the basis of a survey made in late 1913. Prior to this time 
the only official capacity data were published in 1898, 1901, 1904, and 1907 in the directories 
of the American Iron and Steel Association, the predecessor of the Iron and Steel Institute. 
Since the 1914 date capacity figures have been reported annually for steel ingot production, and 
Since 1926 for pig iron production, though capacity estimates for the rolling mills are not pre- 
sented. The reason given for the latter omission is that estimates cannot be made for rolling 








14(Continued) 

Again, the presence of surplus capacity may be stressed by those already in a business in order 
to deter newcomers from breaking into it. 

"Even where the returns are treated as confidential by the collecting agency, many chances 
oferror arise. Trade statisticians who have developed excellent records of supply and demand 
in their respective industries advise us that they have found the measurement of capacity 
impossible to reduce to a generally satisfactory basis. This is partly because interest in the 
subject is relatively new and partly because of difficulties inherent in the problem. It is not all 
uncommon to find that different officers of the same company submit differing estimates of the 
Capacity of a given plant." 

Edwin G. Nourse, op. cit., pp. 23-24. 

I5Cf. J. N. Behrman, “Notes on the Concept of Capacity and the Measurement of Capital 
Intensity," Logistics Papers, George Washington University, 1952, pp. 6-16. 
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mills since the maximum output figure depends upon which of the hundreds of commodities are 
being produced in the mills. No attempt is made to estimate an average or expected bill of 
rolled goods or to use an input item as the basis for capacity measurement; hence, no Capacity 
figures are published for this division of the steel industry. Moreover, when the capacity of 
the steel industry is referred to without qualification, it is the capacity to produce ingots which 
is implied, 16 

On the questionnaire distributed by the Iron and Steel Institute the instructions for 
computing the pig iron capacity of blast furnaces and the ingot capacity of Bessemer, open 
hearth, and electric furnaces are identical. The instructions state that the annual capacity for 
a firm should be computed by multiplying the number of days in a current year by 24 hours by 
the number of furnaces, less the number of hours furnaces will be idle because of relining and 
rebuilding (down time) and holiday shutdowns, times the expected net tons per operating hour. 

Lost time does not include time for fettling, repairs to bottoms, banks, and tap holes, 
hot work repairs, mechanical repairs, and other operating delays. For a shop, the figure com- 
puted before allowances for down time is called gross capacity and the figure after down time 
deductions is called net capacity. !? From our previous discussion it is seen that gross capac- 
ity corresponds to a theoretical capacity, while net capacity is equivalent to a practical capac- 
ity. 

As clear as these instructions may seem, considerable judgment is required to deter- 
mine the expected net tons per operating hour for a furnace, whether it be a blast furnace or 
an ingot furnace. A furnace does not have an unambiguous rated capacity which serves as the 
basis for the expected net tons per operating hour. For example, the maximum amount of pig 
iron which can be produced in a blast furnace during some period of time depends upon the 
quality of the iron ore used and the characteristics of other input materials, especially coke. 
To account for this variability, the rated capacity of a blast furnace is determined not only by 
the physical character of the furnace but also by the character of the available materials. Such 
materials vary greatly in quality. Iron ore in some districts has an iron content of about 60 
percent of natural weight, while in other districts the best ore may have only 35 percent metal- 
lic content. Hence, the rated capacity of blast furnaces is based upon the operating experiences 
of each furnace; two people rating the same furnace may easily derive different figures. More- 
over, a furnace need not be rated at the same figure over time if operating conditions change, 
and clearly two furnaces of the same type with the same physical characteristics need not have 
the same ratings. 

Similarly, for a steel- making furnace the rated capacity is based on its performance 
during a previous full year of continuous operation, and the characteristics of the input again 
are very important. For an open hearth furnace, for example, about 10 to 12 hours are 
required to produce each "heat" or output of steel with a mixture of molten pig iron, scrap 
steel, and limestone. But if some of the scrap is replaced with molten steel from the Bessemer 
process, as is occasionally done, the "heat" is considerably speeded and the rated capacity 
would increase. 

For new furnaces, both blast and steel- making, the initial capacity rating is based on 
past experiences of comparable furnaces operating with similar raw materials and under 





16piscussion taken from "Steel Capacity Statistics Based on Performance," Steel Facts, 
April 1952, pp. 6-7; and Edwin G. Nourse, op. cit., pp. 254-264. 

17American Iron and Steel Institute, Survey of Steel Capacity, Form 202, January 1, 1954. 
For a reproduction of the instructions, see Steel Facts, April 1952, p. 6. 
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similar conditions. These estimates, however, tend to be conservative, and it frequently hap- 
pens that operating experience reveals that a new furnace has been rated too low. !8 

Estimates tend to be conservative even in the rating of existing equipment, particularly 
in regard to the effect of improved raw materials and technological advances on potential out- 
put rates. Heading the list of things that are increasing ingot capacity is coal washing. In 
October 1948, U. S. Steel started a huge washer at Gary, West Virginia, that produces 20 to 23 
thousand tons of washed coal daily. This coal goes to its plants in the Chicago district, where 
the improvement is reported to have added at least 500 thousand tons to that area's steel- 
making potential, but the full extent of this improvement is not reported in the capacity statis- 
tics. Steel officials reportedly feel that such new techniques require study of operating records 
over a relatively long period before their over-all effect can be calculated. In recent years 
Republic Steel has been converting to high-top pressure blast furnace blowing. A company offi- 
cial reported that the new technique added 105 thousand tons to blast furnace capacity in 1947 
and an additional 75 thousand tons to that figure during 1948, but at the same time he admitted 
that these calculations are conservative. Alloy and stainless-steel melting time and quality 
have both been improved by use of oxygen. Reducing time per heat has boosted capacity. For 
the industry as a whole the gain is substantial, but evidently most shops prefer not to talk 
about it. Faster charging of open hearths and electric furnaces is another promising field for 
capacity gains, but apparently little information is available on these improvements. 19 

Some other features of the capacity measures should be noted. The estimates for a 
given year are computed on the basis of the amount of equipment in the previous year plus all 
changes in equipment and improvements which will be made before the coming April of the 
given year. Moreover, all plants which have been idle less than five years are included in the 
capacity figure for the given year. The officials of the Iron and Steel Institute are aware that 
this expedient is not entirely satisfactory, since the number of years a plant has been idle does 
not completely determine its state of repair and potential usefulness. Yet some such arbitrary, 
but consistent, rule is needed to account for idle equipment. Because estimates tend to be con- 
servative and because investment and improvements for the latter three-quarters of a year are 
not accounted for in the capacity estimates for a given year, it is easy to see why the operating 
rate, or utilization ratio, may be over 100 percent. This actually occurred in 1951 when the 
operating rate for ingot production was 100.9 percent.20 

Another factor which may permit the operating ratio to go above 100 percent is the 
allowance for holidays. The present deduction for holidays is four days, even though one or 
more of these days may actually be worked. In the last few years, for example, when steel has 
been in short supply, at least three of the holidays were working days. The rigid allowance for 
holidays thus tends to understate the capacity estimates in some years. 

In addition to capacity estimates, operating ratios for the industry are presented in the 
Statistics. The operating ratio is found by dividing yearly production by the annual capacity for 
the industry. The operating ratio, however, incorporates a lag which was noted previously. In 
the computation of this ratio for a given year, the production figure is that year's total output, 
whereas the capacity estimate depends on the previous year's equipment plus all changes and 
improvements in equipment expected to be completed by April 1st of the given year. In any use 
of the operating ratio it would probably be useful to correct for this bias. 





18Steel Facts, April 1952, pp. 6-7. 
19G. F. Sullivan, "Steel Capacity,'' The Iron Age, January 6, 1949, pp. 198-205. 
20American Iron and Steel Institute, Annual Statistical Report, 1951, p. 7. 
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The capacity statistics show that from 1926, when pig iron capacity was first reported, 
until 1941 the capacity of blast furnaces remained relatively constant, while from 1942 there 
has been a fairly steady increase in reported capacity. The operating ratio, however, reflects 
fairly accurately the course of business activity, and appears to dip each time that the level of 
business falls. The capacity for ingot production, in contrast, increased steadily from 1926 
until 1932, remained constant from 1933 to 1937, and finally increased rather steadily from 
1938 to 1952, except for minor declines in 1946 and 1947. The growth of ingot capacity has 
been relatively greater than pig iron capacity, reflecting the increasing importance of scrap 
steel as part of the charge in Bessemer, open hearth, and electric furnaces. The operating 
ratio for ingot steel production exhibits substantially the same pattern as does that for pig iron 
production, except that the former ratio is rather consistently several points higher, 2! 

The statistics on capacity for the industry represent a complete coverage of all firms 
producing pig iron and steel ingots. We have, however, noted some features of the data which 
indicate that the capacity measures should be treated cautiously and somewhat skeptically. Of 
prime importance is the fact that rated capacities for furnaces depend heavily upon individual 
judgments which are not easily subject to verification. We have noted particularly that there 
is some reluctance to revise capacity figures as innovations occur and operating conditions 
change. According to one author,22 capacity estimates generally tend to be conservative and 
to lag behind improvements in equipment. Other minor deficiencies of the data are that holi- 
days may be worked, even though deductions are allowed for them, and that in the operating 
ratio formula there is a lag between the capacity data and the production figures. Like much 
statistical information, it is not an easy matter to determine what degree of specious or irre- 


sponsible reporting occurs in the steel industry, and I am at a loss to make an estimate of 
this factor. 


2. Petroleum Refining 

For the petroleum industry, statistics on operations are largely compiled by the Bureau 
of Mines and the American Petroleum Institute. The Bureau of Mines publishes its findings on 
capacity estimates yearly in an Information Circular, generally appearing in the latter half of 
the given year, which purports to show the capacity as of January 1 of that year.23 Reports 
have appeared annually from 1918, except for the war years of 1942 through 1945. The Ameri- 
can Petroleum Institute collects some data originally from the refineries, but its estimates of 
refinery capacity are obtained from the Bureau of Mines. The Institute publishes its findings 
in the annual Petroleum Facts and Figures and the weekly Statistical Bulletin. The Oil and 
Gas Journal also presents an annual survey for petroleum refining, and for the years 1942 
through 1945 it is the only original source for capacity data.24 

In the Bureau of Mines annual survey, capacity figures are estimated independently for 
the two major divisions of the refinery: the distillation process and the cracking process. The 
former is reported as refinery capacity and the latter as cracking plant capacity. For each 
category an estimate is made for the plants actually operated during the year, for shutdown 














2lAmerican Iron and Steel Institute, Annual Statistical Report, 1952, pp. 27 and 37. 

22G, F. Sullivan, op. cit., pp. 202-203. 

23For example, cf. F. 5. Lott and J. G. Kirby, "Petroleum Refineries, Including Cracking 
Plants, in the United States, January 1, 1950,"' Information Circular 7578, Bureau of Mines, 
August 1950, 19 pp. 

24For example, cf. R. B. Tuttle, "U.S. Refineries Have Crude-Oil Input Capacity of 6.75 
Million Barrels,'' Oil and Gas Journal, March 23, 1950, pp. 302-304. 
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plants, and for plants in the process of being built during the year. Shutdown capacity includes 
capacity of inoperative portions of operating plants, as well as capacity of plants completely 
shut down. Capacity of facilities under construction includes the capacity of construction to 
augment existing plants, as well as new plants. Before 1938, however, shutdown capacity 
includes only capacity of plants completely shut down, -and capacity of facilities under con- 
struction includes only that of new plants in process of erection. In each circular a time 

series is presented for total refinery plant capacity and cracking plant capacity, according to 
the threefold classification of operating, shutdown, and building. For the given year summaries 
are presented for refineries and cracking plants by districts and by states, and a summary is 
presented for cracking plants by types, i.e., reformed, thermal, and catalytic. A detailed 

report is also offered for each company in the United States where these are classified accord- 
ing to the type of basic atmospheric distillation equipment employed. 

The distillation capacity of a refinery is represented in terms of the input, crude oil. 
According to this concept, the capacity of a plant is the maximum amount of crude oil which 
can be processed in an average day when the equipment is in complete operation. An average 
daily figure is reported, since the capacity figure purportedly shows the maximum daily crude 
runs to stills that can be maintained for an extended period. By this is meant that the daily 
figure allows for necessary shutdown time for routine maintenance, repairs, and breakdowns, 
but assumes no interruptions from strikes, a continuous supply of crude, no material change 
in operating methods, and a demand sufficient to take the entire production. The reporting 
procedure evidently is not as standardized as for the steel industry, since the amount of crude 
oil intake to be deducted for necessary interruptions to continuous refining depends entirely 
upon the individuals estimating plant capacity. Apart from allowances for interruptions, the 
rating of a plant also depends heavily upon individual judgment since the estimate depends upon 
the quality of the crude oil distilled, which may vary considerably, and upon the previous opera- 
ting experience for the plant, which is dependent not only on the size and the quality of the dis- 
tillation equipment but also on the efficiency of the working force and the auxiliary equipment, 
previous price policies, etc. 

The cracking plant capacity is computed in a similar manner and is subject to similar 
difficulties. But whereas distillation capacity is derived in terms of an input, cracking plant 
capacity is computed as the amount of finished product in barrels, primarily gasoline, which 
can be produced in an average day by using charging stocks of the type in current use.25 

For the distillation process, capacity is not measured in terms of output, primarily 
because of a product mix problem. The difficulty arises because of the large number of 
products which are feasibly producible from crude oil by the combination of distillation and 
chemical treatment; and since many products have a seasonal demand, the product mix for a 
refinery varies over the year, even if it specializes in certain commodities. As we noted 
previously, however, the input measure appears to be a reasonably good compromise for these 
difficulties. 

The capacity statistics for the petroleum industry indicate that total distillation capac- 
ity has increased uninterruptedly, except for minor setbacks in 1933 and 1943, from 1918, the 
date of the first annual Bureau of Mines survey, until the present. From 1921, the year oper- 
ating capacity was first reported, it has grown each year, except for minor declines for five 


25The information on capacity, except where otherwise noted, is taken from the Bureau of 
Mines Information Circulars, various issues,and Petroleum Facts and Figures,various issues. 
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years; but, as is to be expected, the rate of growth is less steady than for total installed dis- 
tillation capacity. Until 1943, the capacity of shutdown plants averaged roughly 10 percent of 
total capacity, but with fluctuations between 4 and 14 percent. Since 1943, the capacity of 

shutdown plants has decreased fairly steadily and in 1953 was only 2 percent of total capacity, 

The capacity of cracking plants was first reported in 1925 and until 1938 was computed 
in terms of the charging stock processed by the plants. As nearly as can be determined from 
comparing the two bases of measurement, total cracking plant capacity has grown steadily 
without interruption since 1925, and since 1938 appears to have approximately tripled. During 
this period the relative amount of shutdown capacity has fluctuated from 4 percent to around 
22 percent of total installed capacity. For every year since 1939, however, the capacity of 
shutdown plants has been less than 10 percent of total installed capacity, decreasing to 1.2 
percent in 1953. 

No time series on operating ratios are presented, but the over-all operating ratio for 
distillation plants for the given year is indicated in each Bureau of Mines annual survey. For 
years in which business activity is at a high level the operating ratio between production and 
total installed capacity appears to be around 90 percent. 26 

In the petroleum industry a concept of capacity and interest in measurements of capac- 
ity do not seem to play as central a role as in the steel industry. In the trade journals, e.g., 
Petroleum Processing, Petroleum Refining, and the Oil and Gas Journal, one rarely finds any 
discussion of the measurement of capacity. Generally, any reference made to capacity is to 
report, uncritically, the findings of the Bureau of Mines. In its annual survey the Bureau of 
Mines reports briefly and generally the capacity concept which is the basis of measurement 
for the industry, but it does not discuss either the reporting procedure or the accuracy of the 
capacity estimates. In an interview with officials of a petroleum refining firm, however, the 
author was warned that the capacity measures reported in the various publications must be 
used with caution, It was stated that one cannot place much credence in the accuracy of the 
measures of capacity since the estimates tend to be biased fairly consistently on the optimistic 
side. For many companies it would take major investments for restoring and improving the 
plant before the reported capacity figure could be approached, and apparently some of the 
smaller refineries outrightly and without justification tend to overstate the capacities of their 
plants. 








3. The Electric Utility Industry 

For the electric utility industry, statistics on the production of electric power and on 
the capacity of generating plants in the United States are reported systematically by the Fed- 
eral Power Commission and the Edison Electric Institute.2?7 The data presented by the Federal 
Power Commission is based on required monthly reports to the commission; and the coverage, 
which is substantially 100 percent, comprises around 3700 generating plants operated by some 
1400 electric utilities. Included are plants of both the privately owned electric utilities and the 
publicly owned group. Data on electric energy generated by industrial firms and agricultural 
organizations for their own use are not included in comprehensive form. The capacity of small 
standby plants, operating infrequently and having little or no production, is included in the 





263.G. Kirby, "Petroleum Refineries, Including Cracking Plants, in the United States, 
January 1, 1953,'' Information Circular 7667, Bureau of Mines, September 1953. 

27Federal Power Commission, Electric Power Statistics, and Edison Electric Institute, 
Statistical Bulletin. 
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summary data, though regular monthly reports are not requested from such plants.28 The 
Edison Electric Institute data are derived partly from their own sources and partly from the 
official Federal Power Commission information. 

For the industry various measures of capacity are computed for different purposes. 

The basic data for each measure, however, is the rated capacity, or name-plate rating, of 
generating equipment.29 In brief, the name-plate rating is the manufacturer's estimate of the 
number of kilowatts which can be produced in one hour by the generator, but caution is required 
in accepting the rating as an accurate estimate of the productive ability of the equipment. 

The name-plate rating is the basis of guarantee from the manufacturer to the purchaser 
covering the efficiency of a machine under certain specified conditions. In practice, because of 
differences in operating conditions or of modifications since installation, the actual productive 
ability of generating units may be quite different from the name-plate rating. In their own 
reports most companies carry a "capability rating"' for equipment, instead of the name-plate 
rating, since it more nearly sets forth kilowatt power under actual operating conditions. 

Neither the name-plate nor the capability rating, however, gives the full measure of 
power of generating facilities. Most generating equipment driven by steam turbines or by 
steam engines is capable of considerable overload, often for quite an extended period. With 
generators driven by water wheels or water turbines, stream flow or water flow govern the 
possible output, but even under adverse conditions such generators are capable of carrying a 
substantial overload.30 

With these limitations in mind we may consider the different capacity measures. The 
estimate given as the industry capacity is simply the sum of the name-plate ratings of the 
equipment in service at the end of the year, as reported to the Federal Power Commission. 

To compute annual utilization ratios (commonly called plant factors), an annual capacity is 
obtained for the system by multiplying the name-plate ratings of the equipment in service at 
the end of the preceding year plus one half of the net increase during the year in question by 
8760, the number of hours in a year. This measure thus corresponds to a theoretical capacity 
where potential output is considered as 100-percent use of the installed equipment operating 
continuously throughout the year. 

In comparing peak loads (the highest hourly rate of generation of electric power) with 
capacity, two additional capacity measures are used: dependable capacity and net assured 
capacity. Dependable capacity is defined as the name-plate ratings of all generating plants 
that can be relied upon for active or standby service, plus the capability of hydro plants under 
the most adverse water conditions and plus net firm purchases—all at the time of the annual 
system peak. Net assured capacity is dependable capacity less reported required equipment 
reserves.31 These two concepts, especially the latter, correspond roughly to practical capaci- 
ties. In net assured capacity no explicit allowances are made for maintenance and repair time, 
but the deduction for equipment reserves apparently accounts for this category. An electric 
generator reportedly requires very little repair, but since electric utilities are required to 
meet any demands for power made on the system by users, it is necessary to keep sufficient 





28Federal Power Commission, ''Production of Electric Energy and Capacity of Generating 
Plants 1954,"' Electric Power Statistics, Binder Section 1, 1955, p. 1. 

29Capital equipment in the industry is divided into three major classes: generation, trans - 
mission, and distribution. On the average, for the individual firm there appears to be about 
equal investment in generation and distribution, while investment in transmission accounts for 
somewhat more than half of either of the former categories. 

30Edison Electric Institute, Statistical Bulletin, 1952, pp. 4-5. 

3lFederal Power Commission, op. cit., Report No. 55-66, p. 16. 
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reserve equipment to forestall any power shortage caused by breakdowns of machinery. 
According to the data, reported required reserves presently average around 8 percent of 
installed name-plate capacity. 32 

No information is currently given concerning the policies used to determine required 
reserves, but Nourse maintained that the reserve cannot safely be less than the capacity of 
the largest generating unit or the largest tie-in in the system. He also reported that more 
conservative managers maintain a reserve equal to the largest two units, visualizing the 
chance that some day at the hour of the peak load, when the largest unit is down for repairs, 
the second largest also may suddenly fail.43 The stream flow limitation on potential output 
is becoming less and less important in the United States, primarily because hydro plants are 
continually becoming percentagewise less important in the industry but also because it is 
becoming common to store water power. 

Once the annual capacity is given, the plant factor is computed by dividing annual out- 
put in kilowatts by the capacity. The average annual plant factor for all plants is presently 
around 60 percent, though there exists considerable variation in the system. The main rea- 
sons for this relatively low ratio are easily seen to be the theoretical capacity measure used 
as the divisor and the seasonal demand for electric energy. The capacity concept does not 
allow for the need to withdraw generating equipment from service for periodic overhauling or 
repair, for equipment held in reserve, and for the variation in stream flow for hydro plants. 
As far as the seasonal demand is concerned, the system load varies hourly, as well as daily 
and monthly, but enough equipment must be ready for use to meet peak demands. 

It is interesting historically that from 1920 through 1938 the annual plant factor for the 
system averaged around 36 percent, except for the depression years of 1930-1935. From 1938, 
however, there has been an upward secular trend—the ratio has increased steadily from 35.5 
percent in 1938 to 58.2 in 1953.34 According to the Edison Electric Institute, this improvement 
is the result of several things. The two most important are: (1) the greater reliability of gen- 
erating and transmission equipment which, with integration and interconnection between elec- 
tric systems, has reduced the spare equipment required for emergency or reserve; and (2) the 
constant expansion of uses for electrical devices has prolonged and made more uniform the 
demand for electricity throughout the day. Data on the plant factor extend back to the Census 
of 1902 when the ratio was 25.7 percent and are available at five-year intervals until, beginning 
with 1926, we have the annual record of the Edison Electric Institute. 

The increase in the plant factor since 1902 would appear still larger were it possible to 
allow for changes in the practice of rating electrical equipment. In the early days of the indus- 
try, the typical manufacturer, held to a guarantee of performance and with little experience to 
guide him, seems to have rated his generators more conservatively than today. Users finding 
their machines capable of operating safely at loads greater than name-plate ratings operated 
them at loads above that contemplated by the designers. If plant factors for the early years 
were calculated on the basis of safe capacities, the factors would be somewhat lower than 
reported, 35 

To measure the variations in demand, the electrical industry has developed the concept 
of the annual load factor. This ratio should not be confused with the plant factor, though it is 





32Ibid., p. 6. 

33Edwin G. Nourse, op. cit., pp. 328-331. 

34Federal Power Commission, op. cit., p. xiv. 

35Edwin G. Nourse, op. cit., pp. Sig-321, and pp. 478-480. 
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computed in a similar manner. For a plant the load factor is the total annual kilowatts pro- 
duced, divided by the peak load on the plant (the highest hourly rate of electrical energy gen- 
erated by the plant for the year), times 8760, the number of hours in a year. In other words, 
it is the average hourly production divided by the peak load. For the system the peak load is 
computed as the sum of the peak loads for all the plants, and it is a noncoincident peak since 
the individual peaks may have occurred at different hours, on different days, and in different 
weeks. A true peak, determined simultaneously for the whole system for the same hour, 
would be substantially less than the noncoincident peak.36 The similarity of the noncoincident 
load factor to the plant factor is obvious, the difference being that in the former the divisor is 
the peak load times 8760 while in the latter it is the annual capacity. 

As to the recent accuracy of the data, some evasive reporting may occur or some 
omissions may be present, but the reporting procedures and the characteristics of the equip- 
ment appear to insure that the margin of error is very small, at least as far as the reporting 
of name-plate capacity is concerned. Each piece of generating equipment has a manufacturer's 
name-plate rating in terms of kilowatts, and each plant is required to report these ratings 
monthly to the Federal Power Commission. As noted previously, however, even though we may 
feel confident that the data reported are accurate, caution is required in accepting the name- 
plate capacity as an accurate rating for generating equipment. 


THE USES OF INDUSTRY MEASURES OF CAPACITY 

The uses of industry measures of capacity, of course, are highly conditioned by the 
accuracy of the observations. We have already argued that on this matter the user of capacity 
data must be extremely cautious, particularly since, in many cases, it is very difficult to obtain 
information on the accuracy of the estimates. We need not repeat our previous warnings in this 
regard; instead, we shall confine the discussion to uses of the measures apart from the limita- 
tions imposed by inaccurate observations. 

Capacity measures of this type are, mainly, a function of the existing level of capital 
equipment in physical units and the production relation between services from the capital equip- 
ment and the commodities produced by the industry. Whether or not there is an active interest 
in an industry in actual capacity measures, the information required for the derivation of capac- 
ity estimates is immediately important for the planning of production and investment by indus- 
try entrepreneurs. Other variables which may enter in the determination of capacity, e.g., 
down time and other repair time, also are of importance in this respect. A firm may not 
explicitly combine such data to derive capacity estimates, but implicitly, at least, capacity 
figures of some types are important for the rational functioning of the organization. Conse- 
quently, whether for understanding the formulation of plans or for making predictions, capacity 
measures for industries, where they are obtainable, potentially, at least, embody relevant 
information. 

Generally speaking, it is clear that industry measures of capacity contain useful infor- 
mation. The problem now is what specific uses do they have for economic analyses. Conceiv- 
ably, the estimates may be useful as measures of the ability of industry to produce output during 
the present period of time, as measures of the growth pattern of industries, or as measures of 
capital equipment which may be incorporated into investment models. 





36Edison Electric Institute, Statistical Bulletin, 1952, p. 37. 
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In the first instance, as indicators of the ability to produce for an industry, the meas- 
ures have obvious limitations. The usual industry capacity estimate purports to show the 
productive ability of existing physical equipment, but this ability, of course, is dependent upon 
the full set of conditions occasioning the use of the equipment, and the measure holds only for 
the special conditions which are the basis of measurement. Unfortunately, the majority of 
reported industry capacities are computed on the assumptions that the costs of using the equip- 
ment are not relevant, that there is a free supply of complementary inputs, and that there is 
sufficient demand to purchase any quantity of the commodities produced. These usual assump- 
tions do not realistically portray existing conditions in an industry and are generally not the 
type of conditions which would interest the investigator of the industry's ability to produce. 
For these reasons industry measures of capacity, by themselves, are not very useful for this 
purpose. 37 

As indicators of growth, industry measures of capacity are on a sounder basis if they 
are estimated consistently from year to year or if appropriate corrections can be made for any 
inconsistencies which do appear in the data. They are frequently used for this purpose in trade 
journals and occasionally in studies on the growth pattern of industries.38 In this use they may 
have certain advantages over other measures of growth, such as the value of sales or output, 
balance sheet totals, or the value of durable production equipment, in that they: allow both for 
changes in the amount of durable production equipment in the industry and for technological 
changes, whether these relate to improvements in existing equipment or to the introduction of 
more efficient machinery. Moreover, they may have additional advantages over value meas- 
ures of growth in that they enable the user to show the growth pattern of the industry without 
being troubled by valuation problems. Compared with output time series, their advantage is 
that they are less apt to be influenced by cyclical fluctuations in business activity and con- 
sequently they more readily indicate trend factors which tend to be obscured in the output 
series. On this latter point, Hickman showed, by plotting capacity series, that changes in 
capacity over time seem to be largely independent of cycles in production. The general 
impression seems to be one of secular growth or decline, and this is confirmed when the data 
are studied systematically.>? An obvious reason would appear to be that the firm does not 
vary the quantity of equipment it holds with every short-run fluctuation in demand; this would 
be too costly. instead, it will tend to vary the intensity with which the machinery is used for 
short periods. Another factor would appear to be that technological advances occur even in 
depressed years and tend to raise the capacity rating for existing equipment as well as for new 
equipment. The major drawbacks of the capacity measures, in this use, are that the estimates 
do not indicate the age distribution of machinery and may include an unspecified amount of out- 
put which is attributable to obsolete or long-idle machinery. While capacity measures tend to 
be upgraded for technological advances, it is difficult to account for loss in efficiency of aged 
or obsolete machinery in noncost measures. These limitations may partially account for the 
cyclical insensitivity of capacity estimates. During depressed years firms may continue to 
carry obsolete equipment and may not replace operating equipment as frequently as in pros- 
perous years but, nevertheless, may not downgrade the capacity estimates correspondingly. 





37For a discussion of more appropriate capacity measures for this use, cf. Edward Zabel, 
loc. cit., pp. 237-243. 

38One example is the study by Bert G. Hickman, "Cyclical Variations in Capacity and 
Capacity Utilization in Selected Industries,''*Conference on Research in Income and Wealth, 
National Bureau of Economic Research, Inc., October 9-11, 1953. 

39Bert G. Hickman, op. cit., pp. 12-13. 
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Another drawback in using capacity series as indicators of growth is that it is dangerous to 
make direct comparisons among different industries as to absolute changes in growth on the 
pasis of capacity estimates. The obvious reason is that different industries use different 
concepts of capacity. And this difference in concepts leads to estimates of potential output 
which are. not directly comparable. 

As variables in investment models, both industry capacity estimates and utilization 
ratios may be useful. A capacity estimate may be used as a measure of the existing level of 
capital equipment in an industry, which can be easily related to an output or demand variable 
and which, again, avoids valuation problems. In conjunction with utilization ratios, it is pos- 
sible to set up an investment model which attempts to predict changes in capacity rather than 
directly predict changes in capital equipment. The former type of model would have the advan- 
tages that it avoids some valuation problems relating to capital equipment, the measure of 
capital equipment is more readily related to production, and predicting changes in capacity 
may be simpler than predicting changes in physical equipment in the sense that variations in 
capacity include both technological advances and changes in the amount of equipment. On the 
latter point, for example, firms may meet increases in demand by improving existing equip- 
ment rather than by adding new equipment. Models which predict variations in capacity would 
account for this, whereas the usual model relating directly to capital equipment generally does 
not allow for this type of adjustment. If the ultimate problem, however, is to predict changes 
in capital equipment, the advantage of the capacity model may be nullified in the sense that, 
in relating changes in capacity to equipment, the problems which have been overcome in the 
model arise again in the final transformation. 

Utilization ratios may be helpful in predicting the amount by which capacity increases 
in response to an increase in demand but, as with the capacity estimates, these ratios should 
be used with care. The user must be cautious about making inferences from these data, par- 
ticularly in regard to the absolute levels of the ratios. For example, if industry A on the 
average has a utilization ratio of 60 percent and industry B 95 percent, one cannot infer that 
equipment is being used less efficiently or that more idle or obsolete equipment exists in 
industry B. The difference may be entirely due to different capacity concepts, the accuracy of 
the data, and different operating characteristics for the two industries. Industry A may use a 
"theoretical" capacity concept as the basis of measurement and may include all idle equipment, 
whereas industry B may use a "practical" capacity concept and include only equipment which is 
being currently used or which has been idle for less than, say, three years. In industry A the 
capacity rating of equipment may depend on judgments which would vary among individual 
raters, but in industry B the ratings may rely on verifiable engineering data. Also, industry A 
may be subject to seasonal demand, while in industry B demand may be fairly constant over the 
year, , 

Even for the same industry the user must be cautioned when making inferences from 
year-by-year comparisons of utilization ratios. In discussing capacity estimates, we have 
already noted that over the years the reporting procedures may change in an industry or the 
extent of bias in the estimates may change (e.g., for the electric power industry); these anoma- 
lies, of course, will affect the utilization ratios and should be guarded against. We have also 
noted that in the case of electric power there has been an upward trend in the utilization ratio 
because of improvements in equipment, the practice of interconnection of power systems, and 
some leveling out of seasonal demand. When comparing utilization ratios for the same indus- 
try, particularly over long periods, the user must be alert for such trends. In comparison with 
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previous years, a relatively low utilization ratio may also reflect only the fact that the industry 
is increasing capacity in anticipation of future demand, a healthy sign, rather than validating an 
inference that equipment is being less efficiently employed than in previous years or that the 
industry is relatively depressed. 

In summary, the discussion in this study has revealed the range of concepts which are 
the bases for estimates of industry capacity and has indicated some possible uses for the data, 
To date, however, industry measures of capacity have found relatively little employment in 


economic analyses. 
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The existing body of economic theory on resource allocation rests rather critically on 
the assumption of perfect divisibility of factors. Whenever the factors involved are indivisible 
and are to be combined in small numbers, the principles of marginal analysis offer no guid- 
ance, and entirely new considerations seem to be called for. To obtain an idea of the peculiar 
difficulties that may then be faced, we propose to study the following simple but perhaps not 
unrepresentative problem. 

Suppose that a certain number of passengers is to be transported on a given route for 
which planes of two types are available. Type 1 planes have larger capacity, Si, but smaller 
operating cost, C,, per seat. Provided the total operating cost of the type 2 planes is smaller, 


(1) S9Co < 81C), 


the cheapest transportation for a certain number of passengers may call for the use of one or 
more of the less-efficient type 2 planes. In this paper we obtain a bound on the number of type 
2 planes that should ever be used. This bound is of interest in itself for the question of the 
maximal feasible investment in small planes. But its main usefulness is in reducing the num- 
ber of trial ana error solutions that need be examined in determining the optimal combination 
of planes for transporting a given number of passengers. Computational problems proper, 
however, will not be considered here. 

Suppose the optimal solution calls for the use of n type 2 planes, with possibly some 
number of type 1 planes. It then follows that NS» is not divisible by S}- This is so because 
if NS» = ks, (k an integer), then the passengers in the n type 2 planes could be transported 
more economically in k type 1 planes. 





*Manuscript received April 27, 1956. 
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(2) 


where m is an integer and 1< r< S;- Ba 

We know that the number of passengers in the n type 2 planes can not be transported 
more economically with the use of one or more type 1 planes. In particular, there would be no 
saving in using (m+ 1) type 1 planes. Hence 


NS9Co <(m + 1) SC, . 


Then by use of (2) we have 


—-=t+1]s,c 
S} Jere 


and finally, 


S;- 1 cy 
(3) n< 





So "Cy - Cy 


Similarly, there would be no saving in using m type 1 planes to transport ms, passen- 
gers and using type 2 planes to transport the remaining passengers. Since the number of 


remaining passengers is at most (nS - ms,), the number of type 2 planes needed is not greater 
than 


NS»o - my - Mp * 
53 82 





It follows that 





NS» - ms, So -1 
NSo Co < MS, Cy + 3 + 5 SoCo» 


< MS, Cy + NSpCo - MS) Co + (So - 1) Co - 





INEFFICIENT INDIVISIBLE UNITS 
We then get 


(Sp - 1) Co 


ms, < ——————_ 
1 Co - Cy 


By use of (2) we have 


(Sp - 1) Co 
Co ni Cy 


NSy-1r < ’ 


So - 1 Co r 





n< . + 
Sp fg - Cy Sq’ 





< 


; + 
me “ry % 





(4) n < + 
&g- % S 


Combining (3) and (4), we have 





S,-1l c 
(5) n< Min( - : 


For brevity, write 


Then (5) becomes 
(6) n < Min ( a : +a). 
= 1 ri p > 


From (1) we have a lower bound on p, namely 


(7) 
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1 1 
The first term of the right hand side of (6) is the smaller one when p < a* dy’ and the second 


term is the smaller one for larger values of p. The second term increases with p and 


1 
becomes (s, - 1) when p=1- (Sq - 1) . Since (sy - 1) is an obvious bound for n, (6) should 


not be used for larger values of p. Also, since n is an integer, we can omit any fractional 
part of the number obtained for its bound. Letting [x] stand for the largest integer not exceed- 
ing x, we can summarize the above by 


The bounds (8) increase with both p (the ratio of costs) and A (the approximate ratio of sizes), 
as must be expected for the true maximum of n. 

As an example, suppose the capacities of DC-6's and DC-3's on first class flights are 
58 and 36 passengers, respectively. Then 


19 So i i 1 653 


hema. t= 


12’ s, 29” 2*2.* 38° 1X (Sq - 1) > 665° 


The bounds on the number of DC-3's which need be considered for possible use, for various 
cost ratios, are given below. 
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The problem of finding the optimal number of type 1 and type 2 planes to transport a 
given number of passengers can be formulated as a discrete linear programming problem. 
Let 


Xy number of type 1 planes, 
Xo = number of type 2 planes, 


p = number of passengers to be transported. 





INEFFICIENT INDIVISIBLE UNITS 
Then, we seek integral values of x, and Xq which minimize 
(9) Tr.” wwe 
subject to the constraints, 
(10) X181 + XoS 2 P 
(11) ,20- 


Denote the optimal solution by x? ‘ x). Comparing its cost with that of using 


0 
x X;, + 
7% 3 


Xo Oo, 





(12) 


one can obtain inequality (3) with < for n. 


0 
Xo S 
a 


Also, letting 





one can get inequality (4). 
The linear programming formulation 
permits the following geometric interpreta- 
tion of the allocation problem. In Figure 1, 
combinations of numbers of the two types of 
planes are represented by points of a quadratic 
lattice having integral coordinates. All com- 
binations having a capacity for at least p pas- 
Sengers are located on and above a line, to be MINIMUM COST 
called an isocapacity, whose slope is - s,/s, is ene 
and whose intercept on the x, axis is p/s;. 
The family of isocosts (lines of constant cost) 
consists of lines of greater slope, - C18) /CoSo; ISOCAPACITY 
than the isocapacities. The optimum com- 
bination lies on that isocost nearest to the 
origin which passes through a lattice point 
to the right of or on the isocapacity line for 
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p passengers. This means that the interior of the triangle, formed by the isocost on the right, 
the isocapacity on the left, and the Xy axis at the bottom, should not contain any lattice points, 

If there are more than two unit sizes of the indivisible factor, the inequalities (8) may 
still be applied to yield a bound on the input of any unit in terms of the sizes and costs of larger 
units. Let the units be labelled in order of decreasing size and (hence) of increasing average 
costs, 


then (8) becomes 


(14) n, <Min Min ij 
a: s <j 1 - Pi; 


The constraint i< j is essential, for we must have Pis < a 

Every cost minimization problem suggests an associated maximization problem of 
similar structure. (This should not be confused with the duality principle of linear program- 
ming). For example, in the airplane problem one can ask for a bound on the number of type 2 
planes that should ever be used to transport the maximum number of passengers having avail- 
able a fixed budget. This companion piece will be treated by specific reference to a more 
practical revenue maximization problem. A firm with a limited capital budget considers the 
purchase of two types of equipment having unit costs of q) and oP) and yielding returns of r, 
and To per dollar, respectively, where 9; > 49 and ry >To: (Here qi and Tj correspond 
respectively to S;C; and 1/c i in the airplane problem). Then the best combination that the 
firm can purchase depends on the exact amount of capital available. But we may ask again for 
the largest number of the less-efficient production units that should ever be bought. 

Let n be the maximal number of type 2 production units that should ever be bought. To 
determine n, we need consider only the case that the budget limit q equals nq,. Comparing 


the revenue from n production units of type 2 with that from the use of only type 1 units, we 
have 





INEFFICIENT INDIVISIBLE UNITS 


> = 
ma2%2-\q, |%171° 


nds a, -1 
mas"2-\q, a, / “171? 


NGoro 2 NGor, - (a, -1)r,, 


(15) < 
os G2 Ty - Tp” 





which is analogous to (3). 
Next compare the revenue from n production units of type 2 with that from a combina- 
tion of m units of type 1 and as many units of type 2 as the remaining budget permits. We 


have 


Ndo - mq, 
“i «. i ee 








E - =| ; Ndo- MG, dp - 1 
0) ~ «Ag Q2 


it follows that 
Mq4Ty + Molo - MqyT - (do - 1) To < Mol , 


and then 
Wo - 1 To 


16 ; 
(16) m < aq", 





Now let m take on its largest value compatible with a given q = Ndo. Then 
ndg < (m+ 1)q,-1, 


and from (16) we get 





Fp 2 ce 1 


To q4,-1 
maa ad 
1 2 2 
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Finally, 








+ 
3-3 @ 





(17) n < 








which corresponds to (4). As in the airplane 
problem, here also there is an obvious upper 
be bound of (a; - 1) for n. 
\, The geometric interpretation of the 
solution is similar to the one in the airplane 
\. problem. In Figure 2 all combinations which 
LINE 




















MAXIMUM REVE? 
BINAT! 





cost less than the fixed budget are to the left 
\ 2 aa _—_ of the budget line. We let isorevenue lines 
denote combinations of constant revenue. 
Then the optimal combination lies on that 
isorevenue line which passes through a lat- 
tice point to the left of or on the budget line and leaves no lattice point in the interior of the 
triangle bordered by the isorevenue line on the left, the budget line on the right, and the Xx} 
axis at the bottom. 

If there are more than two types of units, we let 
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Figure 2 


q,-1 
ar i?" iy = 7%. 


and a bound on the number of units of the j-th type that should ever be purchased is given by 


Vij 1 
1 . < Mi sls to a8 ¢ :- ’ 
(18) ny af min( [5 ~ + ‘iz a +44 q; ) 


which corresponds to (14). 
One important difference between the results in these two problems is that in the maxi- 


mization problem the bound given by (15) should be ignored because it is never less than the 
one in (17). Hence (18) can be replaced by 


(19) n. < Min min(| 
as j 





VARIANTS OF THE HUNGARIAN METHOD FOR ASSIGNMENT PROBLEMS! 


H. W. Kuhn 
Bryn Mawr College 





The author presents a geometrical model which illuminates variants 
of the Hungarian method for the solution of the assignment problem. 











1, INTRODUCTION 

The Hungarian Method [1] is an algorithm for solving assignment problems that is 
based on the work of D. Kénig and J. Egervary. In one possible interpretation, an assignment 
problem asks for the best assignment of a set of persons to a set of jobs, where the feasible 
assignments are ranked by the total scores or ratings of the workers in the jobs to which they 
are assigned, It is a special case of the transportation problem, which asks for the allocation 
of a homogeneous commodity from given supplies at a set of sources to satisfy prescribed 
demands at a set of destinations, while minimizing the total source-destination transportation 
costs, The extension of the Hungarian Method to transportation problems can be found in the 
work of L. R. Ford, Jr., and D. R. Fulkerson [2] and of J. Munkres [6]. Comprehensive reviews 
of these problems may be found in the papers of M. M. Flood [3] and T. S. Motzkin [4]. 

The original purpose of this note was to present a modification of the Hungarian Method, 
inspired by an algorithm of Marshall Hall, Jr. [5], for the choice of a system of distinct repre- 
sentatives for a family of sets. But is now seems that a synthesis of all of the available vari- 
ants would be more valuable. Accordingly, the paper now presents a geometrical model that 
is adequate to define and compare the algorithms contained in [1], [2], [5], and [6]. In Section 
2, the Hungarian Method is reviewed in outline so as to isolate the procedure that is altered in 
these variants. The modification suggested by Hall's work is presented in Section 3. A trans- 
lation into graph-theoretical terms of the problem solved by all four versions is the subject of 
Section 4; this section is self-contained and can be read by itself. A final section compares 
these alternative solutions by means of this geometric model. 

2, RESTATEMENT OF THE HUNGARIAN METHOD 
As considered in this paper, the assignment problem asks: Given an n-by-n matrix 


A= (a;,) of non-negative integers, find the permutation iy, Sie i. of the integers 1,...,n 





1The work reported by this note was supported by the Office of Naval Research Logistics 
Project, Department of Mathematics, Princeton University. 
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that minimizes the sum ‘ay; es anj . (If a maximum permutation sum is demanded, set 
1 n 


a= max, a4; and replace A by A' = (aj;) with aij =a- ai, for all i and j.) The Hungarian 
method is then based on two assertions: 

A. The problem is unchanged if the matrix A is replaced by A' = (aj), with aj 2 aij és 
u; - Vj for constants u; and Vip one 1,3= 1, ..., 8. 

B. (Kénig's theorem) The minimum number of lines (rows and/or columns) needed to 
contain all of the zeros of A is equal to the maximum number of zeros that can be chosen, with 
no two on the same line. 

As presented in [1], the algorithm divides into two routines. Routine I is essentially an 
iterative procedure for determining the minimum number of lines and the maximum number of 
zeros referred to in Kiénig's theorem; it is the possible variations in this routine that form the 
main subject of this paper. The input of the routine is a subset of the zeros of A, marked by 
asterisks, with no two asterisks in the same line. If the number of asterisks is n, then the 
desired minimum is zero, and the problem is solved by iy, ee in where the asterisk in rowi 
appears in column ij, for i=1,...,n. Ifthe problem is not solved (i.e., if there are k<n 
asterisks) then the output of a single application of Routine I is one of the following disjoint 
alternatives: ‘ 

Ia. A new set of k + 1 asterisks on zeros in A, with no two on the same line; or 

Ib. A set C of k< n lines which contain all of the zeros of A. 

Routine I is repeated until there are n asterisks (and the problem is solved) or until 
there is an occurrence of Alternative Ib. In the latter event, let h > 0 be the minimum entry 
not appearing in the lines of C. Routine II calls for the addition of h to all entries in each 
line of C (if an entry appears in two lines, then 2h is added) then for the subtraction of h from 
every entry in A. Routine II leaves A non-negative but decreases the sum of all entries by an 
amount nh(n - k). The only property of Ib needed for the finite termination of the combined 
algorithm is that the set C contains less than n lines. The original Routine I produces the 
exact minimum k; in the modification of Routine I proposed in the next section, the set C con- 
tains exactly n - 1 lines at every occurrence of Alternative lb. 


3. MODIFIED ROUTINE I 

Search each column of A in turn for a 0*. Ifa 0* is found, proceed to the next column 
(if no columns remain then there are n asterisks and the problem is solved). If a 0* is not 
found in the column, then the column is called pivotal and is searched for all of its 0's. If no 0 
is found in the pivotal column, then the remaining n - 1 columns contain all of the zeros of A 
(Modified Alternative Ib). If there are 0's in the pivotal column, then their rows are searched 
in turn for a 0*. If such a row contains no 0*, then the 0 in that row and in the pivotal column 
is marked with an asterisk (Alternative Ia). If each 0 in the pivotal column has a O* in its 
row, then these rows are listed in any order: iy, o dtety i,- We now construct a sequence of 
rows based on this initial segment. At step s of the construction, we search row i, for a 0*. 
If a 0* is found, then the row indices of all 0's in its column are added to the sequence in any 
order and any indices already present are omitted. If no 0* is found in row i. then a transfer 
(see [1] and Section 4 below) is possible. Namely, i, was added to the sequence via a 0 inthe 
column of a 0* in an earlier row i, of the sequence. If the asterisk is transferred from row 
a. to row i, in this column, then we have the same situation earlier in the sequence. Con- 
tinuing, we ultimately free an asterisk in a row of the initial segment, iy, eee i, But then 
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an asterisk can be marked in this row of the pivotal column (Alternative Ia). In the remaining 
case, a finite sequence iy, sea iy is constructed with the following properties: Each row in 
the sequence contains a 0* in some column. These columns, together with the pivotal column 
(y +1 columns in all), contain 0’s (marked or not) only in rows iypeees i,. Hence, the 
remaining n - v - 1 columns combined with the rows ij,.-., i, constitute n - 1 lines that 
contain all zeros (Modified Alternative Ib). 

In a rough attempt to evaluate the advantages of this variant, it may simplify the search 
for a transfer in Routine I at the possible expense of smaller changes in the sum of the entries 
in A on the next application of Routine II. 


4, A GRAPH-THEORETICAL EQUIVALENT OF ROUTINE I 

In this section, the problem solved by Routine I (and by the algorithm of Section 3) is 
given an independent statement and then rephrased as a problem in graph theory. 

Let A be an n-by-n square array in which certain of the places are occupied by zeros. 
An assignment is a subset of k zeros, distinguished by asterisks, such that no two asterisks 
lie in the same line (row or column) of A. An assignment is complete if every zero lies in the 
line of some asterisk. A cover is a set of lines of A that contain all of the zeros in A; ifa 
complete assignment is available, only the lines in which asterisks appear will be used in a 
cover. Kénig's theorem asserts that the largest number of asterisks in an assignment is equal 
to the smallest number of lines in a cover. Such an assignment is called maximal; such a cover 
is called minimal. A transfer is possible when there is a sequence of asterisks at (i, j;) \ 

- (i,, i,) in A such that there are zeros at (ig, j,), ee (i,, ip-1) and a zero at (ij, ig)» 

for some io» with no asterisk in its column. The transfer removes the asterisks from 
lip, iy), -- +» (i,, j,) and assigns them to (i, jg), ... » (ip, jp_y)- If there is a zero at 
(i, +p jp), for some i,,,, with no asterisk in its row, then the result of the transfer is an 
incomplete assignment (and another asterisk can be placed at (i r+D i.) ). 

To each complete assignment for A, we shall associate an oriented graph G by the fol- 
lowing rules: 

The nodes V of G are ina 1-1 correspondence with the asterisks. A directed edge 
joins vi to Vo if there is a zero in intersection of the column of the asterisk-node Vv; and 
the row of the asterisk-node V,. The node V is distinguished as a source (sink) if there is a 
zero in the row (column) of V with no asterisk in its column (row). Note that a node can be 
both a source and a sink. 

The equivalence of various concepts in the matrix A and the graph G is displayed in 
the table below: 








Matrix A Graph G 





Transfer Directed path (possibly void) originating at a source. 





Transfer yielding incomplete | Directed path (possibly void) from a source to a sink. 
assignment 





Cover Assignment of + and/or - to each node so that every 
source (sink) is marked +(-), and no edge starts at a+ 
and ends ata-. 
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The first two equivalences follow directly from the definitions. To prove the last 
equivalence, suppose that a cover is given for the matrix A in which the assignment is com- 
plete. Each asterisk is contained in a line of the cover; if it is covered by a row (column), 
mark the corresponding node by a + (-). An asterisk-node V is a source (sink) if, and only if, 
there is a zero in the row (column) of V without an asterisk in its column (row). This zero is 
covered; hence the node is marked with a + (-) if it is a source (sink). Each edge in G cor- 
responds to an unmarked zero, and leads from an asterisk-node in its column to an asterisk- 
node in its row. For this zero to be covered, either the former is marked by a - or the latter 
by a+ (Recall that when the available assignment is complete, only lines through asterisks 
appear in a cover.) To show the other half of the equivalence, consider the set of the rows of 
the asterisk-nodes marked with a + and of the columns of the asterisk-nodes marked with a -; 
under the conditions listed, it is asserted that this set is a cover. Clearly, only zeros without 
asterisks need be considered. If there is an asterisk in the row (column) of such a zero, but 
no asterisk in its column (row), then this asterisk-node is a source (sink) and is marked by a 
+ (-). Hence the row (column) of the zero without an asterisk is in the set. If there is an 
asterisk in both the row and column of an unmarked zero, then there is an edge from the latter 
to the former. Since either the asterisk-node in the row is marked with a + or the asterisk- 
node in the column is marked with a - , the unmarked zero is covered. 


THEOREM 1: Let A be a matrix with a complete assignment and let 
G be the associated oriented graph. Suppose each node of G is assigned + 
or - (and not both), so that every source (sink) is marked + (-) and no edge 
starts at a+ and ends ata-. Then the corresponding cover is minimal. 


PROOF: If there are k asterisks assigned in A, no cover can contain less than k 
lines, The cover corresponding to the given distribution of one + or - to each node of G 
contains exactly k lines. Q.E.D. 


We now prove a purely graph-theoretic theorem that will form the basis of our discus- 
sion of Routine I (and, incidentally, prove Kénig's theorem). 


THEOREM 2: Let G be a directed graph in which two disjoint 
(and possibly void) sets of nodes have been distinguished as sources and 
sinks, respectively. Then exactly one of the following alternatives holds: 

(a) There is a directed path in G from a source to a sink. 

(b) Each node of G can be marked with one sign (+ or - ), so that 
every source (sink) is marked + (-) and no edge starts at a + and ends ata -. 


PROOF: Both alternatives cannot hold since, for any distribution of signs and any 


directed path from a source (+) to a sink (- ), there must be an edge that starts at a + and 
ends ata -. 


Let S be the set of nodes that are connected to some source by a directed path (starting 
at the source). If a sink lies in S then (a) holds. If a sink does not lie in S, then assign + to 
the nodes in S and - to the nodes not in S. This marks the sources and sinks as required by 
(b). Suppose that an edge starts at a node marked +, and hence in S. If the directed path from 
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a source to this node is extended by the edge in question, the terminus is connected to the 
source by the extended path. Hence the terminus is in S and is marked +. Q.E.D. 
Remark 1: A proof of Kiénig's theorem now results as follows: It is a simple matter 
to construct a complete assignment of asterisks to the zeros in A (say, by assigning an aster- 
isk to the first zero in each row which is not in the column of an asterisk that has been placed 
previously). If alternative (a) of Theorem 2 holds, a transfer yields an incomplete assignment 
and another asterisk can be placed. Repeating, alternative (a) can hold only a finite number of 
times (surely no more than n - k). When alternative (b) first occurs, we have a minimal cover 
with the same number of lines as there are assigned asterisks, and Kénig's theorem is proved. 
Remark 2: If there are no sources (sinks) then the theorem is trivially satisfied by 
marking all nodes - (+). 
Remark 3: The construction of the proof can be altered in the following manner. Let 
§ be as above. Let T be the set of nodes that are connected to a sink by a directed path 
(ending at the sink). If there is a node both in S and T, then (a) holds. Otherwise, mark the 
nodes in S (T) with +(-), and mark all nodes that are in neither S nor T with the same sign 
(either all + or all -). 


5. VARIANTS OF ROUTINE I ON THE GRAPH G 

The methods described below all amount to a construction of the set S or T in different 
manners. They will be defined in informal geometric terms, and the translation back into terms 
of the matrix A will be left for the reader. 

1. (Kuhn [1]). In the search for a path from a source to a sink, start at any source and 
continue as far as possible along a directed path. If a sink is not reached, backtrack to the 
nearest branch point and try again. When the directed paths originating at one source are 
exhausted, try the next source. If all sources are exhausted and no sink has been reached, the 
set of all nodes that have been encountered is S. 

2. (Ford-Fulkerson [2]). Start simultaneously from all sources and construct all 
directed paths with one edge (or less) originating at the sources. Extend these to all directed 
paths with two edges (or less). Continue fanning out from all sources until a sink is encoun- 
tered or until S is exhausted. 

3. (Hali[5]). Pick a column of A without an asterisk and call it pivotal. If there is no 
0 in this column, then the remaining columns form a cover of n - 1 lines. Otherwise, each 0 
is in the row of a source. Fan out from these sources (as in 2, above) until a sink is encoun- 
tered or until the set of nodes connected to these sources by directed paths is exhausted. The 
cover then consists of the rows of the asterisk-nodes thus encountered completed by all col- 
umns except the pivotal column and the columns of asterisk-nodes already covered by rows 
(n- 1 lines in all). j 

4, (Munkres [6]). Mark all nodes both + and -. Then examine the nodes that have 
two signs, according to some preassigned order (say, from left to right in the matrix A). If 
such a node is a source (sink) or is connected from (to) a node marked only + (-) by a directed 
edge, remove the - (+) from the node and prime the edge used in the criterion (if any). If both 
Signs are removed, a unique path from a source to a sink consisting of primed edges has been 
found. Otherwise, the examination is repeated until no more removals are indicated. At this 
point, the nodes marked + (- ) constitute S (T) and the remaining double signs can be changed 
to -. (See Remark 3, above.) 
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Through the use of the graphical model, the reader can easily concoct variants of his 
own. Moreover, simple graphical examples show that none of the four methods given above is 
the best for all matrices A. 

The method of Munkres sketched in Section 5 above is replaced in the revised version of 
his paper 6 by the following variant: 

5. (Munkres 6). Mark all nodes - . Examine each source, sink, and edge in some pre- 
assigned order (since each of these objects is associated with one or more unassigned zeros in 
A, the natural order of these zeros can be followed). If a source marked - , or an edge connect- 
ing a + to a - , is found then change the - to a + (and prime the edge). If a sink marked + is 
found then a unique path consisting of primed edges from a source to a sink has been discovered, 
(We now note that, if the rows and columns associated with the current marks + and - are 
covered, each of the objects for which we are searching is associated with one or more uncovered 
zeros in A, and we can follow the natural order of these zeros.) Otherwise the process is 
repeated until all sources are marked + and no edge connects a + to a - (Alternative (b) of 
Theorem 2), 
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The fundamental problem of system analysis is to decide which is 
the best system of a given collection of systems. A precise formulation 
of this problem is made and a method of solution given which is, toa 
large extent, independent of the criterion employed. 











1, INTRODUCTION 

Suppose there are given several different systems, each of which serves the same gen- 
eral purpose. With respect to this purpose it is assumed that there can always be associated 
with each system a measure, V, of the value received when this system is employed, and a 
corresponding measure, C, of the cost. How does one decide which of these systems is, in 
some sense, superior to the others? This problem customarily has been resolved by the fol- 
lowing procedure: A criterion D is selected, which is a function of V and C that appears to 
be particularly appropriate to the specific class of systems involved; and it is assumed that 
system i is superior to system j, providing D(i) > D(j). For example, the specific criteria 
most frequently found in the literature are: 

1. D= V for (a)C =Co, or for (b)C Cp; 


2. D=-C for (a) V = Vo, or for (b) V >Vo; 


3. D=V-C; 


4. D=-Vv/c. 
All systems of the class are then compared under D, and that system having the largest value 
of D is called optimal (i.e., it is in this sense that one of the systems is considered superior 
to the others). 

Certain aspects of the above approach appear unsatisfactory, namely: (1) to the authors 
knowledge, no significant correspondence between classes of systems and criteria has been 
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established, (2) there is no possibility of arriving at an optimal solution composed of a mixture 
of the systems from the given class, and (3) it is frequently true that the numerical values of 
several of the parameters that enter into D are entirely unknown. There may also be other 
objections. 

This paper attempts a characterization of the classes of systems for which maximi- 
zation of D results in a function (preference function) that decides which systems are superior 
to the others independent of the specific manner in which D depends on V and C. Incidental 
to this characterization will be the result that for many systems the preference function obtained 
depends on substantially fewer parameters than appear in D. 

To the best of the authors’ knowledge, the first paper concerned with such a characteri- 
zation is by R. W. Shephard (to be published). For the particular case of weapon systems 
appropriate to air defense (a generalization of example 6.1.2, below) and for the specific cri- 
teria 1(a), 2(a), and (4), Shephard has found a preference function by the use of methods dif- 
ferent from those employed here. 

We now summarize the results of the paper for the case in which a decision is to be 
made between two different systems and similar results are true for n different systems, 
This summarization is not mathematically precise, but it will aid in understanding the treat- 
ment which follows. We begin with an example. , 

Assume there are two types of weapons that may be used against a target of initial 
value Vo: Suppose that each strike with a weapon of type 1 (or type 2) is composed of a cell 
containing ky (or kK») planes, only one of which carries a weapon. Also, suppose that p (or q) 
is the probability that a weapon of type 1 (or type 2) will not destroy the target and that Cy (or 
Co) is the cost of a weapon of type 1 (or type 2) plus the expected cost of delivery. These num- 
bers will, in general, depend upon the cell size. The problem is to determine under what con- 
ditions one type of weapon is superior to the other. As pointed out previously, one usually 
compares these two systems in the following manner: carry out one strike with each weapon 
type, measure the value V and cost C associated with this strike, choose a criterion D(V,C), 
and if D, > Dy then type 1 is superior to type 2. On the other hand, the procedure followed in 
this paper is the following: carry out x strikes using weapons of type 1 and y strikes using 
weapons of type 2, calculate the value V(x, y) and cost C(x, y) associated with these strikes, 
and choose any criterion D. By treating x and y as continuous variables (some of the impli- 
cations of this assumption are given in 6.5), we use elementary calculus to determine the 
behavior of the surface D in terms of the functions V(x, y), C(x, y). In fact, if Vy V., C,, 

C. >0,* it is shown in (3.1) that along either of the curves, V(x, y) = Vo or C(x, y) = Co, the 
behavior of the surface D, for every function D with Dy >0 and Dc < 0, is determined by the 
sign of the function V,/C,, -V /Cy. If this function is positive, the larger value of D will 
occur in the direction of increasing x; if it is negative, the larger value of D will occur in the 
direction of increasing y. Consequently, one can use this function to determine the set of pairs 
(x, y) on which the maximum value of D must occur. If the pairs (x, y) are restricted to a 
particular set R, then the maximum value of D for every function D will occur either on the 
boundary of R or on the set of pairs (x, y) for which V,/C, - V,/Cy = 0. In the following, the 


pair of functions (V,/C,, V,/Cy) is called a preference function. 














*If f= f(x), ...,%,), then f, is the partial derivative of f with respect to x. . 
1 n x; P i 
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Returning to the example, we have 


V(x, y) = Vo(1 - p*a¥) 


C(x, y)= C,x+Coy 


and V,/C, = Vop'a” |1n pi/C,, Vy/Cy = Vop*a lin q\/Co, and if Vp), C,,C, >0, then the 
sign of the function V,/C, - V,/Xy is determined by the function |1n p|/C, - |Inq|/Cy, which 


is independent of D, x, y and Vo- Consequently, if the only restriction in x and y is that they 
are positive, and if |1n p\/C, > |in ql/Co, then the maximum of D will occur for y = 0 and 
some x. In this case, we say that weapon type 1 is superior to weapon type 2. For this 
example, V,/Cy > V,/ey for every x, y, and, in such a case, the function (V,/C,, V,/Cy) is 


called an invariant preference function. General classes of value and cost functions which lead 
to invariant preference functions are discussed in Section 4; in particular, see Theorem 4.6. 


2, FORMAL STATEMENT OF THE PROBLEM 

For the purpose of this paper, a system §$ is defined whenever a rule is established by 
means of which it is always possible to decide whether or not a given object belongs to the 
system. The collection of objects that make up a system fall into natural categories S; ‘ 

8, ... , according to their function or their physical identity, or both. This separation into 
categories is accomplished by using various rules, so that it follows that 8; ‘ So goes gare 

also systems; these are called subsystems. Each subsystem can frequently be separated ina 
similar fashion so that one arrives at a hierarchy of subsystems. The fundamental problem of 
system analysis is to select at every stage of a hierarchy that particular collection of subsys- 
tems which results in the best (optimal) system in some sense, 

For simplicity, the discussion will be started by specializing to a two-stage hierarchy 
of subsystems. Let each subsystem S;; i= 1, 2,...,N, be characterized by the parameters 
s; = (5° baie ; ) so that S, = S; (s;), i =1, 2,...,N. Construct the system § = § (S,, iene 
Sy. Xyo-++ 9 Xn T) from x; identical subsystems of kind S;, i= 1, 2,...,N, where r = (rj, 


wey a, represents the collection of all parameters which affect $ but are exterior to all the 
subsystems S;. Further, it is assumed that there is associated with every N-tuple x = (x, " 


Ms Xy) a measure of the value V(x, S, r) and a measure of the cost C(x, S, Fr), where S = 


(s;, kas Sy). Finally, let D(V, C) be an arbitrary function, increasing in V and decreasing 
in C. We shall introduce the following definitions: 


2.1 Definition. Let U,,---,U, be n variables. 
1. R" is the region defined by u; Es 


2. ¢"™ is the class of all functions defined and continuous over R". 


3. Ci is the class of all functions H(u,,... , u,) such that H, Hy. ae oo 


4. C5 is the class of all functions H(u,,..., u,) such that He cy and Hy, >@ i= i, 3, 


“ 
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C3 is the class of all functions H(u, ee u,) such that He ch and the equation 
Se * u;,0,..., 0) = Hp has a solution uw’, i-=1, 2,...,n for every 


Ho€[0, + wo). 


. ¢* is the class of all functions H(u, » Up) such that HEC : ‘ Au, >0O and *., <0. 


2.2 Definition. The system Xp is superior to the system xX (in symbols Xo >x) if! D(Xp)> 
D(X); %p is not inferior to X (in symbols Xp 2X) if D(X) > D(z). 


Let the functions d, (x, 8), i=1, 2,..., N, be defined for a known set of parameters 
i 


5 and all xeR. Further, let a class {D} of criteria be given, and suppose R'cRN 





2.3 Definition. The N-tuple (4s,> or dg.) is a preference function with respect to R' and 


{D}, if for every Xe R' at which dg # dg for at least one pair i, j, there exists a point 
i j 


%, €R' such that X, >X for every De{D}. 





2.4 Definition. The N-tuple ds er %, is an invariant preference function with respect 


to R' and {D} if it is a preference function for which there exists a permutation of 
(1, 2,... , N) such that for all XeR’, dg dg, >... ds. 


Within the framework of the preceding discussion, the theorems which follow charac- 
terize a class of systems {$} for which the preference function (or invariant preference func- 
tion) is independent of the specific functional dependence of D on V and C. It should be noted 
that V and C need not be measures of value and cost in the following. The discussion is 
merely a treatment of the behavior of a certain class of surfaces {D[V(x), C(x)]}. 


3. PREFERENCE FUNCTIONS 

With respect to the theorems to be proved, the components of X will always be treated 
as continuous variables. By means of an example (Section 6.5), some of the difficulties encoun- 
tered if the components of X are considered integral will be illustrated. 


3.1 Theorem 
2 n 
If* V,Ce Ca, then the N-tuple (Vx, /x,? Vx,/exy» Poe Vi! Cx) 


is a preference function with respect to R' and {D|De ¢*}, where R’ is either of 
the following: 

1. Any subset R' of any member PF of the family of surfaces { V(x) = 
Vo} U {C(x) = Co}; where R' is open in rf (Interior of RN), 





lin the following, the symbol D(X) represents D(V(x,3,F), C(X,5,T)), and where convenient 
the dependence of V and C on § and fF is omitted. Further, the distinction between the parame- 
ters 5 and F is not essential, and they will all be considered as represented by the single 
symbol &§. 

2a€B if the element a belongs to the set B. ae AUB if either a€A or a€B. a€ ANB if a€A 
and a€B. AB if each element of the set A belongs to the set B. 





DECISION PROBLEM FOR A CLASS OF SYSTEMS 


2. Any set R' which can be represented as a union of sets defined in 1, 
and in particular any open subset of the interior of RN. 


PROOF: Let V, Ce c%, and suppose IP is one of the surfaces {C(x) = Co} (for r e { V(X) = 
Vo} the proof is the same). It is sufficient to let Xp be any point of R'c I for which Vx/Cx, > 
v, ie, for at least one pair i, j. Since Cec J the function “x; 7 0 at the point Xo: It then 
i 


J 
follows from the implicit function theorem that the relation C(x) - Co = 0 uniquely defines x; 


as a continuous function of Xipeeey Xj_p x; tpces > *y throughout a certain region in RN 
about Xo and, consequently, throughout a set R" containing Xo and open in R'. Furthermore, 


ax,/2 Xi = a everywhere in R". Hence, for every De c*, 


= Dy [Yx, Vx; (-Cx,/ c,.)| = Dy _ [Vx,/ Cx, ~ Ya" c;,] 


i 


everywhere in R", and D. (x,) >0. Since D is continuous in X throughout RN and R" is open 
x; 0 


in R', a point K€ R" exists such that x? X for every Dec*. Since Xo is an arbitrary point 


of R', (Yx,/xy enna ¥,,,/© is a preference function with respect to {D|De¢*} and R' as 


*) 


defined in hypothesis 1. 
If R' is a set as defined in hypothesis 2, the theorem is clearly true, because through 


any point in R' there passes a surface of the family {C(x) = Co}. 


3.2 Corollary 
If R' is a subset of the 2<m<N dimensional set (*, >0, j = 1, 2, 
J 


o% ba then (Vx /C, ,...-,V, /C,.\ isa preference function with 
x x x 
iy iy iin im 


respect to R' and {D|De ¢*}, where yoy eee XE (Ky, - + 5 My) 
1 1m 


3.3 By inspecting the proof of Theorem 3.1, it is clear that it need only be required that 
Dy >0, De <0 [or Do <0, Dy 2 0], provided that © is restricted to the family of surfaces 


{C(x) = Co} [or {V(X) = Vo}]. This permits the discussion of functions D which are independent 
of C [or V]. 


3.4 Further, P need not have been restricted to curves of the form V = Vo or C = Co- It is 


easily shown that the result is the same when YP is any curve having a slope always lying between 
the slopes of these constant value and cost curves. In fact, this is the largest class of curves Tf 
for which the theorem is true. Since the additional variety in the sets R' which can be admitted 
is not greatly increased, use of only the constant-value and constant-cost curves will be con- 
tinued. 
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The preference function obtained in the preceding theorem can be used to reduce the 
possible sets on which the maximum of D can occur. In the following theorem (3.5), R' is any 
set satisfying property 1 or 2 of Theorem 3.1, i.e., R' = Up (R'N'r) where I ranges over an 
appropriate subset of {V = V)} U {C = Cp}. The discussion is also restricted to two variables 
x and y, and it is assumed that every curve [I is intersected by the curves V,/c,, - Vc, = 0 
at most a finite number of times. 


3.5 Theorem 
For every D(V, C)€C* with V, Cec2, every maximum3 of D on C1(R’) 
occurs at a point (x, y), which lies on one of the curves V,/c, = V,/<y or on 
the set¢ Cl{Up [Cl (T'NR') - (TNR')}}. 


PROOF: Since the curves B(x, y) = V,/c, - V,/Cy = 0 intersect f in at most a finite 
number of points, it follows that these curves divide the set R'NT into a finite number of sets 
Yq» Such that each y, is open in f and the function B has fixed sign (but may vanish) on y ~ 

Continuity of the functions B and D onT implies their continuity on Cl (q) for each a, 
and since B has fixed sign on 7,4, it has fixed sign on Cl (v4). Consequently, the maximum of 
D on Cl(y,,) must lie on the set Cl (y a) - Yq» which consists of the endpoints of Cl (y q) Since 
the set of a's is finite, it follows that the maximum of D on Cl(R'NT) must lie on the set 
Uy, [Clty a -/ ‘al =Sp. It is clear that these points will be either zeros of B or boundary 
points of R'. 

Furthermore, R' = Up (rf R') implies that for any (x, y)¢R' there is an (Xp; Yo)» depend- 
ing on (x, y) and contained in UpSp, such that D(x, Yo) > D(x, y). But this implies that the abso- 
lute maximum of D on Cl(R') must belong to the closed set Cl{UpSp}. Therefore, since each 
set Sp can be written Sp = Cl(rMR') - rNR'+TPu (B = 0), it follows that Cl {UpSp} = 
Cl {Up [C1(P NR’) - rNR'}}+ Cl {Up [PN (B = 0)]}. The observation that C1{Up [IN (B = 0)]} 
lies on the curves V,/c, - V,/Cy = 0 concludes the proof. 


3.6 Theorem 3.5 can be sharpened somewhat. It is clear from the proof that the maximum of 
D can occur only on certain of the curves B(x, y) = 0, but not on all of them. It is only neces- 
sary to consider those portions of these curves where B changes sign along curves at the 
intersection points. Even this class is too large and can be reduced further by elimination of 
the places where a change in sign of B implies a minimum rather than a maximum. 


3.7 To illustrate more clearly the implications of the above statements, suppose the curve 
V,/c,, > Vy/Cy intersects each of the curves C(x, y) = Co in, at most, one point, and suppose — 
these curves are similar to those shown in Figures la and lb. If V,/C, <(>) V,/C, for x< (>) 
X45 then the maximum of D in R? must occur on the hatched lines in Figure la. If V,/C, >«) 
V,/cy for x < (>) X,, then the maximum of D in R2 must occur onthe hatched lines in Figure lb. 


3.8 Theorem 3.5 may be generalized to the case of N variables. In such a case, the conclusion 

is that the maximum of D occurs for some X which belongs either to the set Cl {Up [C1(rR’)- 

(rN R')]} or to one of the sets of points where either dg = dg . dg = dg a ds Waicscig 
.=_ 4 ¥ 


dg = dg, = «+= dg 


3In unusual situations, there may be more than one maximum. 
4By C1(A) is meant the closure of A in R2. In most cases, the set mentioned inthe theorem 
will consist merely of those points which are added to the boundary of R' by closing R'. 
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Figure la Figure lb 


4, INVARIANT PREFERENCE FUNCTIONS 

In this section, the invariant preference functions (see Section 2.4) for a two-fold hier- 
archy will be treated, and for simplicity the discussion will commence with V and C defined 
only over the two-dimensional plane of variables x, y. 


4.1 Even though, in general, the preference function depends on (x, y), it may happen that 
Vc, - V_/C,, does not change sign for any (x,y) in R“. In such a case, it follows from 


Theorem Pd 1 that the pair (V,/C,, v,/c ) is an invariant preference function, and it is clear 


from Theorem 3.5 that the maximum of 5 lies on the y-axis or on the x-axis, according as 

V/C,S Vv /Cy or V,/C, 2v,/Cy, respectively, for all (x,y) in R. Thus, with no side con- 
ditions on x and y, an invariant preference function leads to a preference, independent of the 
functional dependence of D upon V and C, of a pure system over all other systems. This is 


an important class of systems and therefore warrants further study. 

If V,Cec?, the pair (V,/C,,V. /C..) is an invariant preference function if, and only if, 
the expression ViCy - VyC, = f(x, y) #0 for all (x,y). To describe a class of functions V,C 
which lead to invariant preference functions, it would be sufficient to solve this underdeter- 
mined partial differential equation for some f(x, y) #0. This means C could be considered 
an arbitrary function of x,y, and V = V(C,x,y) determined in such a way that the above equa- 
tion is satisfied. Classes of solutions can easily be obtained in this manner, but given two 
functions V and C, it is usually easier to show f(x, y) #0 throughout R2 than to verify the 
required functional dependence of V upon C, x, and y. 

However, some classes for which it is easy to verify this functional dependence may be 
obtained by finding those functions V,C which satisfy the relationship V,/c, = (b/a(V,/Cy), 
Le, av Cc - bv, C, = 0. This differential equation is equivalent to the statement that the 
Jacobian of the functions V(x, y) and C(bx, ay) vanishes throughout R2. Consequently, for every 
given C, a function W exists such that V and C are functionally related through w[V(x, y), 
C(bx, ay)] = 0. Conversely, it is clear that V = U[C(bx, ay)] satisfies this differential equation 
for arbitrary functions U. There are two subclasses of these solutions which seem to be of 
interest, namely, the solutions which correspond to functions V and C which satisfy either 
(1) V/V =b, c,/c = a, or (2) V,/c, =b, v,/c =a. It is easy to show that V = W(bx+ y), 
C = U(ax'+ y), with W and U arbitrary, are the general solutions of (1); and V = bg(x) + ah(y), 
C = g(x) + h(y), with g and h arbitrary, are the general solutions of (2). In the remainder of 
this section, (1) will be generalized to the case of N variables. 
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4.2 Definition. The functions V and C belong to the class PN ona set R'SR if v,cech 
and, whenever their partial derivatives are written for all X€R' as 





Vx, > {4 G7) F, 5, 7) and 


Cx, ns 8; (s, T) G; (x, 8, T), 


where the functions f; and g; are all positive, the integers i= 1, 2,...,N can be partitioned 
into 1< k <N distinct nonempty subsets I, such that for n= 1, 2,..., k: 
1. Fi ~ F; and G; = G; for all (i, j)e ‘ and 


2. fy (ay Big(n) >1,/é; for some ip(ne 1, and all i # ig(n)e 


The proof of the main theorem on functions V,Ce pN 


mas, the first of which may be found elsewhere.5 


requires the following two lem- 


4.3 Lemma 
A necessary and sufficient condition that n functions of n+ p independent 
variables be connected by a relation which does not involve these variables is 
that every one of the Jacobians of these n functions, with respect to any n of 
the independent variables, should vanish identically. 


For the following lemma it is convenient to introduce the notation o., = » £,X; ‘ 
ie 
n 


vn" pe g:x;,n=1,2,...,kand¢= 6,,...,6,),V* (Wy)... ,V,), where I, is defined 
n 


in Definition 4.2. 


4.4 Lemma 
Suppose V(x, soe% Xn) € cN and each ¢,, n = 1, 2,..., kK contains at 


least two of the variables X4,+++» %y- Then the Jacobian of the k + 1 functions 


1, ace 4, V with respect to any k + 1 independent variables % chad ‘ , 
+ 
is either identically zero or of the form A (‘ Vv, -£,V, ); where A isa 
i j 


constant and the indices i, je(i,, er ig 41) nt, for some n. A Similar result 


holds for the functions y ores Vy, and some function C ech! : 


PROOF: It is clear that 26 ,/2x; = 0 if i, #1, and yn, - \ if i,el,. Suppose 
there is a permutation of the indices iy, ---» ip, Such that (‘y ww dg ta) el, (ta, ep 


ta») Ely, .-+ 5 (Ny, _ ~oners k+1) €I,. If some I, contains none of the indices ij, then 





5Goursat, E. J. B., A Course in Mathematical Analysis, New York, Ginn and Co., 1904, 
Vol. I, p. 56. 
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choose @, = @,_4 and omit I, from the above sequence. Furthermore, abbreviate ij by 
j,j=1,2,...,k+1. Then 


fy eee ty, 


0 ($4, «+ 2%: V) 


5 aR. AE 
8 (Ky, +--+ » Xe 4) ; 
1 f aT fet 


-»V 








wi 
i 


where all blank spaces in this determinant represent zeros. Expand this determinant along the 
last row and consider the cofactor of Ve . This cofactor, 
1 


. a 


a,’ 
(- 1)K+ 1 


fog tec? Ake 








is zero except when @, = 2, A» = 3,..., % _ ,=k. If these conditions are satisfied, then 


8 (by, +++ 5 dys V) 
Se ee : fee 


ray Vx,’ Vx,’ ee hs ‘ 





k+1 
=(- 1) 7. ¢,.- F Brews t 
(Vx, {1 Xo 2) f3 k+1. 








This discussion for Vv, holds true for any Vv, , and the desired result is obtained. The proof 
i 


is similar for ¥,,...,¥, and C. 


4.5 Theorem 
If v, Ce PN ona set R'CRN (see 4.2), then there exist functions U, W 
such that U(V, ¢) = 0 and W(C, wy) =0 where the functions ¢, y are the same as 
those in lemma 4.4. Conversely, if V = U(¢) and C = W(W) where U, We ch but 
are otherwise arbitrary, then V,Ce pN on 


PROOF: Let ay set I, mentioned in Definition 4.2 contain at least two variables. If 
the functions V,Ce pN on a set R’ cRN then they must satisfy the system of equations 
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fV, - fV, =0 
) x; ix; 

i,jeI,,n=1,2,...,k. 
gC, - gC 


}~x; xj 


Conversely, if V,Ce ch and these equations are satisfied, then V,Ce pN. for, if v, = fe (x), 
i 


x; j°x; 


Vx. = fn), fy; >0, then f,V, - f,V, =f)f)(£-n) #0 unless Ve pN similarly, Ce PN. 
j j 
The remainder of the proof for the case where each 4 contains at least two variables 
follows directly from Lemmas 4.3 and 4.4. If some of the sets I, contain only one variable, 
the theorem is proved in the same manner by applying Lemmas 4.3 and 4.4 only to those sets 
‘. which contain at least two variables. 


If Definition 4.2 is combined with Corollary 3.2, the following theorem is immediate. 


4.6 Theorem 
If V,Ce PN on RN and De C*, then the tuple ({f,/g;}, ie I) is an 
invariant preference function over any set R'<{XeR |x; fixed for all i¢ 1} 
which satisfies the hypotheses 1 and 2 of Theorem 3.1. 


4.7 If V,Ce pN (see 4.2) on RN. and De C*, then there is a system x*e RN having all compo- 
nents zero, except x (n) [with i,(n) as defined in 4.2], n = 1, 2,..., k, such that X*>x for all 
oO 


xe RN, 
In this manner the N dimensional problem is reduced to a k <N dimensional problem, 
and the preference function in this k dimensional space is * 


fg, 5 wes 
“ig()  “ig(1) 


v, <«& 
*io(n) X; (0) 


4.8 If n = 1, then the n-tuple (f;/g;5 re fy/En) is an invariant preference function over RN, 
If, in addition, f,(s,r) = f(s;,T), g,(S,r) = g(s;,r), i= 1, 2,... , N, then the subsystem S;(s;) 
which maximizes the ratio {(s;, F)/g(s;, F) with respect to s; is not inferior to any other sub- 


system. Consequently, the system design may be based on this ratio. ~ 
I 
5. SEVERAL STAGES 
The foregoing discussion covered exclusively the problem of a two-stage hierarchy. 
This section will indicate how a three-stage hierarchy can be treated by a brief glance over 
the case where a particular invariant preference function exists in the second stage of the 
system. 


Suppose (see Section 4.8) that (4s,>4s,) = (fs /8s fs, /es,)> where fg. = £(5,;T) and 
gs. = g(S;;F), i= 1, 2. Since *, i. implies S, superior to Sp (in the sense that an x, 
exists such that D(x, 0) >D(x,,Xp) for all (x1, xo)€ R?), maximization of the ratio fs /€s, 


will determine an optimal system for the two-stage hierarchy. Hence, to determine the optimal 
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combination of subsystems for a three-stage hierarchy, it is sufficient to determine the com- 
bination of subsystems of S, that will maximize f, /g, and then determine the number of 
1 s, S; 


subsystems S, that will maximize D. Suppose then that 8, = S,(s,) consists of y, identical 

subsystems of kind T, = T,(t,) and of Yo identical subsystems of kind To = T (ty). Then 

5, =5,(¥y,¥qrty,ty, m), where m is the collection of parameters that affect S, but are 

exterior to T, and Ty, and the functions f, and es, depend upon the variables y,,Yo and 
1 H 


the fixed parameters t,,t),m, andr. The problem now is to maximize the function *%, = 
5 (¥yYq)/ Bs, (y,, Yo) with respect to Vy and Yo: This situation is exactly analogous to the 


preceding problem with fs, and gs, taking the roles of V and C, respectively; thus, if it is 


assumed that fy Bg € 4 the same kind of results can be expected. 
a» 2 


It should be noticed that the actual quantities of the optimal subsystems can be deter- 
mined for all stages of the hierarchy except the second, since the complete criteria for these 
stages are precisely known and are independent of D. 


6, EXAMPLES 
6.1 V = V(p*q’), C = C(C,x + Coy). 


If V = V(p"q¥), V' <0, 0 <p, q <1, C =C(C\x + Coy), Cy,Cy >0, C’ >0, then V,, Vy» 


C,,C >0, and, furthermore, V = V(e* Inp+y in 4) = U(x|in p|+y{lnq|). Therefore, it 


y 
follows, from Section 4.8, that the pair 


(1) (jin p\/C,, [In q\/Cy) 


is an invariant preference function with respect to R2 and {D|De C*}. It should be emphasized 
that to obtain the preference function it is sufficient to know the form of the variables on which 
V and C depend (e.g., it is sufficient to know that C depends on x and y only through the form 
Cix + Coy; it is not necessary to know that C equals C,x + Coy), and knowledge of the actual 
functions is entirely unnecessary. Further, it is interesting to notice that the preference 
depends only on p, q, C,, and Cy, independent of all other parameters which may enter into 

the expressions for V and C. In the examples which follow, the complete functions V and C 
can easily be written. In very complicated cases, however, this is not always practical, even 
though it is sometimes easy to reason merely that V and C are functions of p*q’ and C,x + 
Coy, respectively. 


6.1.1 Offensive Bombing Campaign 

This example was discussed in the introduction. We remark here only that it is now a 
Special case of 6.1. It should be noticed that the system which is superior to all other systems 
could be obtained by maximizing the ratio |1n pl/c, with respect to the parameters determining 
Pp and c,. For example, if p and Cy depend only upon k, and the characteristic parameters 
of the weapon, one could determine the optimal cell size for a given weapon by maximizing the 
ratio |In p|/c, with respect to k,. 
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6.1.2 Air Defense 

Assume that the defense of an installation of initial value Yo is carried out by employ- 
ing x weapons of type 1 and y weapons of type 2 against each attacking aircraft in such a way 
that each weapon affects only the aircraft for which it is intended. Assume further that the 
attack is made in K strikes, each composed of a cell containing k aircraft, only one of which 
carries a weapon. The problem again is to determine the circumstances under which one type 
of weapon is not inferior to the other. The value V expected to be received from the expendi- 
ture of x weapons of type 1 and y weapons of type 2 is defined to be the value of the instal- 
lation expected to be saved, plus the value of the planes expected to be destroyed. If p (or q) 
is the probability that a weapon of type 1 (or type 2) will not destroy an enemy aircraft, and if 
Cc, (or Co) is the expected cost of a weapon of type 1 (or type 2), then V = V(p*q") and C = 
C(C,x + Coy). 

Therefore, the preference function is the same as in the preceding example, but the 
expressions for p, q, C,, and Cy will be different. 


6.2 V= V(p*q’), C = C(u* v) 
If V = V(p*q’), C = C(u* vY), v' <0, C'< 0, 0<p,q,u, v <1, then Vxr VyrCy,C, > 0 


and, furthermore, V = U(x|1n p| + ylin q|), C = W(x|1n up| + y/1n v|). Therefore, it follows 
from Section 4.8 that the pair (1n p/In#, 1n q/1n v) is an invariant preference function with 
respect to R* and {DiDe C*}. 


6.2.1 Antisubmarine Campaign 

An aircraft carrier is used as a home base for dispatching cells of x +y +1 planes in 
search of enemy submarines. One of these planes carries search equipment, x of the planes 
each carry a single weapon of type 1, and y of the planes each carry a single weapon of type 2. 
Each cell goes its own way and searches for submarines independently of the other cells. The 
problem is to determine under what conditions one type of weapon is superior to the other. 

The object is to destroy submarines for the purpose of saving the shipping and shore 
installations which the submarine would attack if it survived. In carrying out this program 
there is involved a certain cost, namely, the carrier may be lost because of accidents of planes 
carrying weapons on take-off or landing, and the carrier and land storage base are both sub- 
ject to destruction from predetonation of weapons in storage. It is now desired to determine 
the specific manner in which V and C depend upon x and y. 

The value received from destroying submarines is the difference between the value 
V(x, y) of ships, shore installations, etc. (all targets which submarines can attack) expected to 
remain if the planes attack submarines with x and y weapons on each mission and the value 
V(0, 0) expected to remain if no weapons are employed against submarines (this is independent 
of x and y). The expression V(x, y) is determined by the number of submarines, the number 
of times N any particular submarine is encountered during a cruise, the probability of survival 
A(x, y) during each encounter, and the value V, the submarine is expected to destroy between 
encounters. The function V(x, y) is therefore a function of A(x,y), since N depends (in a com- 
plicated manner) upon x,y only through the function A(x, y), and V,, is independent of x, y. 
Furthermore, the probability of survival A(x, y) during each cilialaunen is given by A(x, y) =P xq) 
where p (or q) is the probability that a weapon of type 1 (or type 2) does not destroy the sub- 
marine. Finally, it is seen that V = V(x, y) - V(0,0) = U(p*q”) = W(x|1n p| + y|1nq|), where W 
is (in principle) a known function. 
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Now consider the cost C. The expected carrier cost is the value of the carrier and its 
contents, multiplied by 1 minus the probability, Py that the carrier survives. If Pa is the 
probability that the weapons do not predetonate in storage and Pi (x, y), P(x, y) are the proba- 


bilities of successful take-off and landing, respectively, then P. = PyP,P fad where m is the 


number of take-offs and n is the number of landings. Assume that n= m6, where 6 is the 
probability of not detecting a submarine on a mission. Further, Pa is assumed independent of 
x and y, on the basis that in either of the storage areas the total number of weapons does not 
vary much during a cruise and the variation in predetonation probability from weapon to weapon 
is small. If, in addition, Hy (or v1) is the probability that a weapon of type 1 (or type 2) is 
loaded on the plane and takes off without a detonation, U5 (or Vo) is the probability that a plane 
carrying a weapon of type 1 (or type 2) lands and is unloaded without a detonation; then Pp = 


6 6 
ui vi, P. = ue vy, and the expected cost of the carrier is a function of (Hy Mo) (v, vo). If 


the expected cost of the land storage base is assumed constant, and the cost of planes lost by 
predetonation while on a —— is ignored » comparison to carrier cost, then the total 
expected cost C = C(x|1n Hy Ho| +ylin Vy v5), where C is (in principle) a known function. 


Finally, from ipetion 6.2, the . hee preference function for this problem is given by 
{in p/1n (Hy 9), 1n g/in (v, v9)I- 

In the preceding example, it is of particular interest to observe that the weapon choice 
is independent of those parameters which would be almost impossible to evaluate. For example, 
the cost of carrier and contents (including lives) and the value of the destruction a submarine 
could accomplish if it survived do not appear in the preference function. 


6.3 V= v(E, ja ax, - Pk *),c -o( E c.% 


If v' >0,C' >0,a,,C, >0, 0<p, <1,i=1, 2,...,N, then V,,C, >0,i=1, 2,...,N. 
i? “i i Xi? “Xj 


Therefore, it follows from Section 3.1 that the N-tuple 


N x 
(2) (4s, ds,» eee 4s.) = & |1n p,| p. P; Vie 7. x| 


is a preference function with respect to {D|De C*} and any set R' satisfying the conditions of 

Theorem 3.1. Furthermore, when R'= Interior (RN) one can conclude (see Section 3.8) that 

the maximum of D must occur for some X which belongs either to the boundary of RN or to the 

set S of all points where either d, = ord, = =d, ,or,...,0rdg =do =...= dg. 
Sj “s si %s; it Sy 


The relation d, = dg implies 


N X, N X, 
a; + lIn pl Py, C;= a; + |in P;| Rath P, C; . 


N 
1 Xe |1n py| = - In [(a;/C; - a,/C,)/(|1n p;|/C; - [An pi|/C;)} . 





JACK K, HALE AND RONALD L,. REED 





Figure 2a 


According to Eq. 3, the conditions dg = 
i 
4g, (i, j= 1, 2,...N, i # j) define a family 


of parallel hyperplanes. From this it fol- 
lows that dg = dg = ds i #j+#k is never 
i j 








satisfied, so that the maximum of D must 
occur either on one of the planes defined by 
(0,0,x)) Eq. 3 or on the boundary of RN. 

Specific values of the parameters 
entering in this problem will in general 
allow one to restrict the maximum of D to 
a set smaller than that mentioned above. 
To illustrate this, suppose N = 3 and the 
set of points where ds. ds, is not empty, 


and the set where % = 4g, is not empty (therefore, qs, # ds, throughout R3). Assume the 


Figure 2b 


plane C(x, , Xo, X3) = C,, is divided as shown in Figure 2a, where the sets A,, Ag, and Ag do not 
include the intersections of the planes %, = dg, and ds, = dg, with the plane C(x) = C o If 


in A 


~ > . in A, and ~ > 4s, in Ay, then it follows that d 


<d UA, and d, <d 
S, os 8, ~ “8 


S3 
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in A,UA,. Since ds, t ., in C(x) =C,, and ds, > ds, > ds, in Apo, it follows that s, > ds, 


in C(x) =C,. Consequently, the system (0, 0, X39) > (x1, Xo, X3) for every (Xj, Xo, Xa )E A,UAs 
and (X19 0, 0) > (x;, Xos X3) for every (x, Xo, Xa )€ A,. Since this is true for every C,, the 


maximum of D must lie on the hatched lines in Figure 2b. 


6.3.1 One realization of the types of V and C functions discussed above can be found in 
Example 6.1.1 if the costs of the weapons used by the enemy are taken into account. 


6.4 Subsystems in Series 


el Boe 
6.4.1 Consider a system composed of x paralleled - 
qT YW | 


subsystems of kind A in series with y paralleled sub- ae 
systems of kind B, as schematized in Figure 3. 

One realization of such a system would arise Lg (3) 
ina circuit where B represents a timer which, ata 
prescribed instant, allows a switch A to be actuated. 
Asecond realization relates to an air defense problem, 
where the assumption is made that it is possible to stop an enemy attack only if at least one 
radar warning device and at least one antiaircraft weapon function successfully. 

Returning to the general system, suppose A has a failure probability p and a cost c 
and that B has a failure probability q anda cost d. If it is assumed that the value V expected 
to be received is proportional to the probability that the system functions and that the cost c is 
the sum of the component (subsystem) costs, then 








Figure 3 


V=V,(1- p) (1-9), 


“C= xc +yd. 


According to Theorem 3.1, the pair (V,/C,, V,/Cy)= [v,( -q°)p*|in p \/c, V(1 - p*)q) |1n q|/d] 


is a preference function. By Theorem 3.5, the maximum of any De C* lies on the curve V,/Cc, - 


Vy/Cy = 0 or on the boundary (x = 0 or y= 0 in this case), where 


V,/Cy - Vy/Cy = (Vop*a® |1n p|/c) [aY - A(p* - 1) - 1), 


= c|inq|/d|inp|. 


If (x,, y,) is any point on the curve q'Y - (p™* - 1) - 1= 0, then for (x = Xo» ¥ >Yq), V,/C, - 
V,/Cy >0O and for (x >X 5» Y= Yo)» V,/c,, - V,/Cy <0, since this curve passes through (0, 0) 
and has a positive slope. Hence D increases as the curve V,/c, - V,/Cy = 0 is approached 


from either side along a constant-cost or constant-value curve. This means the maximum of 
D lies on the curve V,/C, - V,/ey = 0 and not on the boundary of R?, The character of this 


curve is shown in Figure 4 for the case C, >Cy and |1n p|<|1n ql. 
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6.4.2 It is instructive to go a step further with 
this type of example and allow the designs of A 
and B to vary. Suppose X subsystems of type 
A, and Xo subsystems of type A», all in parallel, 
and in series with a parallel arrangement of Xy 
subsystems of type B, and x, subsystems of type 
Bo (see Figure 5), are introduced into the system. 
In the order thus established, let the failure prob- 
abilities be Py» Po» Iy> Ig and the costs Cy, Co, 
d;, do. Under the same assumptions as made in 
the preceding example, 








x, x Xq X 
V=V,(1- Py) Po”) (1 - a> a9*) , 
C= Xy Cy + Xq Co + Xq dy + Xq do. 
Figure 5 Without loss in generality, it is assumed that |1np,| /e> 
|In pp /eo and |1n 9 |/do > |1n q; 1/4}. It follows that 


V and C belong to the class p4 on the space R‘ of vectors X= (x1,Xo, Xa, X4), since the index 


integers 1, 2, 3, 4 can be partitioned into two sets: (1, 2) = I, and (3, 4) = I,, with F, = Fo= 


VoPy Po (1 - 4; 4 ), F, es F, at Vo G9 (1 - Py Po ), f, = |1n P,|, f,= |1n Pol; f, = |1nq)|; 


{,= }1n Qo |, 81 = Cy, Bo = Cos Bg = Gy, B4 = do, {,/e; > fo/ and f4/&4 > f,/e3. Consequently, 
Theorem 4.6 applies, with the result that (|1n Py!/cy, jin Po /e9) is an invariant preference 
function over the subspace of vectors (x, Xo, 0, 0), and (|1n q,|/d,, |In Qy|/do) is an invariant 


preference function over the subspace of vectors (0, 0, Xa, X 4): This means there always exist 
x, and x, such that (x,, 0, 0, x4) >(x,', X9', X3', x,') for any system (x,', X»', X3', X4'). Thus, 


the system has been reduced to all A,'s and all Bo's, so that, according to Section 4.7 
(Vx /C, » Ve ) is a preference function, and one can proceed as before to obtain a 
+. i. 


relationship between X, and X4, as in Figure 4. 


This example brings out an interesting feature of systems corresponding to V, Ce pN 
namely, optimization can be carried out in each subspace (i.e., corresponding to I) independ- 
ent of the other subspaces. In the above example, this means the designs of A and B are 
indepefdent of each other and can be based on max(|1n p|/c) and max(|1n q|/d), respectively. 


6.5 The Variables x and y Integral 
In general, the assumption that x and y can assume a continuous range of values yields 
a solution which is not an integer pair. For many applications (in particular, those included in 


this paper) such a solution is physically inadmissible. Furthermore, the solution may not cor- 
respond to the nearest integer pair. 
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To illustrate the restrictions imposed on the various parameters by requiring x,y to 
be an integer pair, consider (in connection with Example 6.1.1) the particular function 
D(V, C) = C subject to the condition V 2 Vo i.e., minimize C = (C,x + Coy) subject to 


v= vo(1 - p*q’) 2V> or, therefore, subject to x|In p|+ y|Inq| >A, where A = [1n(1 -V,/Vo)I- 


Geometrically (Figure 6), this means the line C,x + Coy = C is to be moved downward and to 
the left, subject to the requirement that at least one point on this line shall lie above and to the 
right of the line x|1n p| + y|/Inq|= A. 


Region of Admissible Values of (x, y) 


x{In p| + ylin q| = A 





Figure 6 


If x and y were continuous this would be a simple linear programming problem, the 
solution to which is (x,, y,) = (0, A/|1n q|) or (x,, y,) = (A/|1n pl, 0), according as C,/C, 2 


\In p|/|In q|. This result is the same as obtained previously in Example 6.1.1, as would be 
expected. 

However, if the solution (x, y) is required to be an integer pair, the above analysis is 
valid only if y, = |1n(1 - V,/v,)|/l1nq| is an integer, or x, = |1n(1 - V,/v,) |/\1n p| is an 
integer, according as Ci, 2 |1n p|/|1n q|. In this connection it may be possible to adjust Vo 
somewhat to make y o OF X, come out integers, but the implications of such an adjustment are 
not easily assessed. If Xo and Y, are not integers, then the solution may be on an axis or on 
the interior of the admissible region, depending on the relative magnitude of all the parameters 
P, 4; C., Co, Vo and Vo: The integral case will not be considered further. 


6.6 Maxima and Minima of Functions 

Preference functions can sometimes be used advantageously to discuss maxima and 
minima of functions of two or more real variables x, y,... . As a simple example of the 
procedure, consider the function 





2 
D(x, y) = (xX + by + ¢) ;x>0,y>0. 
x" +y +1 


So that the correspondence with the preference process shall be clear, the identifications 
V = (ax+ by + c)? and C =x? + y- +1 are made. Clearly Dj, is negative, and Dy, V,, Vy 
C. C. are all positive; and all functions involved and their first partial derivatives are con- 
tinuous, so that it is permissible to employ the previous results. 
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B(x, y) = V,/c,, - V,/Cy = [(ax + by + c)/xy] (ay - bx) , 


the maximum lies on ay = bx, or x = 0 or y=0. For y >(b/a)x, B >0, so that D is increas- 
ing as the curve ay - bx = 0 is approached from above along, say, a constant C curve. For 

y < (b/a)x, the function B <0, and D increases as ay - bx= 0 is approached from below, 
along the same curve. It follows that if D has an extreme value in R2, it is not a minimum 
and it lies on the curve ay - bx = 0. In this way, the problem is reduced to finding the maxi- 
mum of 





D( Bs) [x (2+ =) 


Xx, =3) = 
at b? 
x 1+-> 
a 


a 


which can easily be found by standard procedure. The result is: D has a maximum of 


a? +b? 


+ c located at the point (x, y) = (a/c, b/c). 
There are two attractive features of this method: it may save needless computation 
by showing the most direct way to arrive at the solution, and it makes the sufficiency easier 


to establish. 


7. CONCLUDING REMARKS 
7.1 After the notions of preference function and invariant preference function were introduced 
in Sections 2.3 and 2.4, it was established that (Yx,/x) ioe Ving! Cx) is a preference 
function. Some of the interesting features of this result are: (1) this function is independent 
of the specific dependence of the criterion D upon V and C; (2) it is unnecessary to know the 
specific functions V and C, it being sufficient in every case merely to know the relations 
among the components of X entering in V and C; and (3) it selects from the sometimes vast 
number of parameters entering in V and C only those essential for the selection of a pre- 
ferred system. This latter feature is particularly important since many of the parameters of 
D may be entirely unknown. 

If V and C belong to the class of functions pN (see Section 4.2), then the preference 
function (v \ ee a ) takes a much simpler form, and very specific results can 

Ky Xy Xn XN 


be obtained (see Section 4.6, et seq.). 

Two concluding remarks (Sections 7.2 and 7.3) concerning the interpretation and use- 
fulness of the above-mentioned results are followed by some comments on open problems 
(Section 7.4). 


7.2 Problems in system analysis frequently come encumbered with side conditions of remark- 
ably diverse character. When this happens it is necessary to inquire into their origin and 
nature for the purpose of determining to what extent they can be relaxed or entirely ignored. 
Such an inquiry is important because the analyst must be certain that he is actually finding a 
best system, in some sense, rather than merely verifying a group of opinions by admitting too 
restrictive side conditions. 
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Nevertheless, it will sometimes occur, even for systems of the class discussed, that 
there are one or more essential side conditions. In this event, the preference function 
a ay, may frequently be used to discuss the behavior of the D surface 
my * *y *N 


inthe region R' of interest. In fact (see Theorem 3.5), it was shown for N = 2 that the maxi- 
mum of D must lie either on the 'boundary' of R' or on one of the curves V,/c,, = V,/Cy- 


1.3 In the case of an invariant preference function over RN) it is always possible to decide if 

a system is not inferior to all others which are to be compared. But no information is obtained 
concerning the number of subsystems which should be employed (i.e., the number that leads to 
a maximum of D), because this number depends specifically on D and thus will be, in general, 
different for every D. If it could be shown that some D, say D,, has greater significance than 
any of the others, the number of subsystems could be determined by using D,. If sucha D, 
exists, then it could be maximized from the outset, obviating the necessity for finding a pref- 
erence function. However, as long as there is uncertainty about the relative merits of the 
various admissible D functions, it may be that the only sense in which one system can be 
preferred to another is in the sense of 'superior' or 'not inferior’ for all criteria D, as defined 
herein. 


1.4 The only fundamentally unrealistic assumption made was that of continuity of the compo- 
nents of X. It would be highly desirable to obtain analogous results when the restriction to 
integral components of X is strictly adhered to. However, even the most simple example (that 
of linear programming) introduces great difficulties. 

Finally, an interesting extension still within the assumption of continuity in x would be 
to consider the vectors X as strategies in a two-person game where the payoff is given by any 
member of a certain class of D functions. For example, suppose the players share a common 
value V(X, ¥) but have separate costs C,() and C,(y). In this case, D can be any well-behaved 
function of V, C,, and Co, which is increasing in V and Cy and decreasing in C,, where D is 
considered the payoff to the player corresponding to X. If a saddle point of D exists, then 
theorems analogous to those proven in the foregoing can be established. 


The authors wish to thank Sandia Corporation for permission to publish this paper, and 
Dr. R. W. Shephard for allowing them to read his paper prior to publication. 








ON THE THEORY OF REPLACEMENT OF MACHINERY 
WITH A RANDOM FAILURE TIME! 


George H. Weiss 
U. S. Naval Ordnance Laboratory, 
White Oak, Silver Spring, Maryland 





This paper studies the reliability function, and moments thereof, 
of a machine which is periodically replaced in order to increase the 
mean time to failure. Two cases are studied: strictly periodic replace- 
ment and randomly periodic replacement. 











IL INTRODUCTION 

Some time ago N. R. Campbell [1] presented an analysis of several replacement poli- 
cies for continually operating elements subject to random failures. His treatment led to 
several difficult integral equations. As a result, it is difficult to draw any qualitative con- 
clusions from his work. 

Specifically, Campbell considered the question of how many replacements are necessary 
to keep a series of street lights in continual operation. He assumed that the failure of each 
light could be detected instantaneously and that instantaneous replacements could be made. 
Under these assumptions Campbell treated two replacement policies, the first replacing the 
lights as they failed and the second replacing all of the lights periodically. 

This paper will consider the problem from the point of view of time to first system 
failure with a given replacement policy. We will consider two possible replacement policies: 
a strictly periodic policy and a randomly periodic policy. 


I. NONRANDOM PERIODIC REPLACEMENTS 

We treat first the case of strictly periodic replacement. Let R(t) be the reliability 
function, i.e., the probability that a device survives for a time t or greater. If R(t) is dif- 
ferentiable, let p(t) = -dR/dt. Let us assume that a given mechanism is replaced every to 
units of time and inquire for the probability density for a failure at time t. 

Let R,(t) be the reliability function if k replacements have been made and no failures 
have occurred. Then, if kt <t<(k+ 1)to, we have 





1 This work was performed while the author was employed at the Ballistic Research Labo- 
ratories, Aberdeen Proving Ground, Aberdeen, Maryland, 


*Manuscript received March 14, 1956. 
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(1) R, (t)= [R(ty) * Rit - kt). 


The n'th moment of this distribution, written <t%, is given by 


e (k+ 1)to 


<t>=e- 5 t" dR, (t) 


(k+ 1)to 
[Rit,)]* J t" dR(t - kto) 


0 0 


t 
0 
k 
- 5, [Ri(ty)] i (x + ktp)" dR(x) 


t 
k 2 in ee ee 
» [Rttoll Ps (j) (to) f x/d R(x). 


The summation over k can be effected in closed form by manipulation of the geometric series. 
For n an integer we have 


x n 
3 of ht = 5 eo St 
(3) i*s (x 4) 1-x 


where the a are Stirling numbers of the second kind [2]. Denote by «ti > o the integral 
s ty j 
(4) <tl>, =- { x) dR(x). 
0 


Then the n'th moment <t™> can be written 


ak k 
~~ + Ae Ww G,_j k! [Rlto)] 
(5) <t">= Pe (7) th <td>» Z kei 





1 [1 - Rito) 


The first moment or expected time to failure is given by 


to R (to) +<t> 
1-R (to) 


t 
[° Rwat 
_ 4 


= 1 - R (to) 





<t>= 
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and the second moment is given by 


to to 
2 { tR(t)dt 2ty R(ty) { R(t) dt 
0 P 0 
[1 - Rity)] [1 - Rit)I? 





<t?>= 


A point of some interest is the limiting form of these moments when to ~0. If we 
suppose that R(t) has at least one derivative at t= 0 then the integrals in Eqs. (6) and (7) can 
be approximated for tp ~0 

to 
t R(t) dt “> 


to 
{ R(t) dt ~ty, 
0 


1 - R(t) ~p (to) to ’ 
so that for to ~0 we have 


1 
<t>~ 
p (to) 


02 acts ~ tp? ~_1 _ 5° 
[pe (to)] 


The limiting form of the first moment indicates that an indefinite increase in the expected time 
of first failure can be achieved by shorter and shorter replacement intervals only if the first 
derivative of the reliability function is zero initially. 

In [3] Herd gave a brief discussion of replacement policy based solely on the criterion 
of failure rate. He considered two cases: a monotone increasing failure rate and a monotone 
decreasing failure rate. When an increasing failure rate describes the mechanism a periodic 
replacement program is advisable, while in the contrary case a periodic replacement policy 
would remove from the system mechanisms with desirable failure characteristics. The value 
of<t>can be used in a similar fashion. If<t>is a decreasing function of to then a replacement 
policy would be beneficial. If<t>is an increasing function of to then no replacement policy is 
advisable. An intermediate case is also possible. Consider the reliability function 


1 2 
(10) R(t)=1- 2 - ar : 
(1 + at) 





; *Herd defines failure rate as ". . . the probability of failure in a finite interval of time, 
given the age of the product." 
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The function<t>is found to be 





2,2 
1 {((1+ a@t,) (3+ 9at, + 5a“tp) 
(11) Wale cone G 0 ei. 
3a a ty (2 + ato) 


A plot of 3a@ <t>as a function of aty is given in Figure 1. It can be seen that<t>has a mini- 
mum which occurs at a@ to ~ 1.29. The curve in Figure 1 has an asymptotic value equal to 5, 
The function<t>is equal to its limiting value for to ~0.33 a" “, so = it is clearly inadvisable 
to institute a replacement program unless the condition to < 0.33a@°~ can be met. 















































Figure 1 


We can also examine this case by using a criterion derived from that of [3]. The prob- 
ability that a given mechanism will fail sometime during the time interval (h, h + t), conditional 
upon its survival till time h, is R(t|h) where 


R (h) - R(t +h) 


(12) R (t}h) = R (h) 





If we define the hazard rate, H(t), as3 er R (dt|t)/dt, we find 
dt> 


(13) Hi) =- am ort, 





3This is the "force of mortality" as defined by Davis in [4]. 
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as given by Herd. For the particular reliability function given in Eq. (10), the function H (t) 
is given by 


4a (2at+ at?) 


14 H(t) = . 
_ (1+ at) (1+ 4at + 2at?) 





Agraph of H(t)/4q@ as a function of at is shown in Figure 2. This function has a maximum 
for t~0.5a@~°. The hazard rate increases until t reaches this value, then decreases there- 
after. Both of these criteria—that in Eq. (11), and Herd's which is essentially Eq. (14)—would 
suggest a replacement program if it is possible to make replacement faster than a certain 
rate. A more comprehensive analysis of the factors governing the choice of a reliability 
program would take into account the economics of the situation. A discussion of replacement 
program from this point of view is to be found in [5]. 






































at 


Figure 2 


A particular calculation of some interest concerns the value of <t>for a network with 
n components in parallel, each with reliability function R(t) = exp(-a@t). The reliability func- 
tion for the network is 


(15) R(t) =1- (, a Rie | é; 


A calculation carried out in the manner of Appendix A of [6] shows that 


-at 
0 
B, (1-e ) 


‘on. 





(16) a<t>= 


l-e 
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where the functions B,, (x) are 


2 3 n 
e a es 
(17) B,(x)=x+5 +3 + += 


-at -at,\" 
For a single component we have a<t>= 1 so that the functions B, (1 -e 0) / (1 -e 0) 


represent the improvement afforded by the use of a replacement program and redundant com- 
ponents, Several curves of a<t>as a function of aty are shown in Figure 3. For large n the 
functions B, (x) can be written 


(18) B, (x) = - log (1 - x) + 0(x"), 
|x| <1 


so that Eq. (16) becomes, approximately, 


to 
(, - . to)" 


Equation (19) indicates that when a policy of finite replacement is adopted, the mean time to 
first system failure goes up approximately exponentially with the number of redundant compo- 
nents, If there is no replacement, or equivalently exp(- ato) = 0, we have 


(19) <t>~ 


(20) B (1) 21444... ¢ 2 -10gn, 
n 2 n 

















By (I-e72to )(1-e"A to? 






































Figure 3 





REPLACEMENT OF MACHINERY 


where the last estimate follows from the well-known asymptotic expansion of the harmonic 
series (cf. [7]). Hence, if no finite replacement policy is adopted, the effect of adding redundant 
components shows up approximately as the logarithm of the number of components. 

We can also examine the distribution of the number of replacements made before a sys- 
tem failure occurs. The probability that a given device (or system) will function correctly for 
atime ty is R(t). Hence, the probability that there will be n replacements before a system 
failure occurs is [1 - R(tp)] [R(ty)]”. If the moments of this distribution are defined by 


(21) <n*>= [1 - R(to)] am nk [Rit))]", 
n= 


then we have 


R (to) 


ade ——— 
1- Rito)’ 


R(t 1+ R(t 
<n’> = ( 0? : es F) 
[1 - Ritg)] 





R (tp) 
o? =<n?>-<n>* = a . 
[1 - R(to)] 


In our previous paragraphs we have allowed the reliability function to depend only upon 
the time. In general, one should study the reliability function as a function of other parameters, 
such as temperature, humidity, or other relevant environmental factors. With this in mind we 
may write the reliability function in the general form, 


(23) R(t)= R(&,, Bo, -.. By ; t) 


where now the 8's may be random variables. Let us now define F(x,, Xo, --- Xj; t) as the 


distribution function of the #'s at time t, that is, 
(24) F(x), Xo, coe * ; t) =Pr{p,< X1»B9 £ Xo, +++ By $x at time t} ° 


In order to get the values of <t"> we may define an averaged reliability function R*(t), where 


@ 
(25) R= ff Ra, d--- 2,5 ta @ a F(a, Oy... 0,5 t), 
ca 1% +++ “n 


and then take moments with respect to R*(t) as in Eq. (2). As an illustration, let us consider a 
reliability function R(t) = exp(-a@ t), where @ has the probability density k[exp(-ka)]. The 
function R* (t) is 
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@ 
(26) R* (t) = k f en Ahk +t) gq= *, 


which implies that for periodic replacement <t>is given by 


(27) <t>=k (+=) In (1 3 


We can see from Eq. (27) that even though the original reliability function is an exponential the 
first moment is a function of to- 


Ill. RANDOM PERIODIC REPLACEMENTS 

It can happen that it is impractical to replace a mechanism in a strictly periodic fash- 
ion. This might, for example, be the case when a given mechanism has a variable work cycle 
and is of such a nature that a replacement in midcycle is impossible or impractical. In this 
eventuality the replacement policy would have to be a random one, taking advantage of any free 
time available for making the replacements. , 

To describe the replacement policy we will assume that there exists a function F (t) 
which is the probability that no replacement is made inthe time interval (0, t). If it is certain 
that some replacement will eventually be made, F(t) satisfies 


(28) F (0)= 1, F(w)=0. 


A further condition on F(t) is that it be a monotone nonincreasing function. It is also possible 
to describe a periodic replacement policy in terms of F(t) by taking 


F (t)= 1, t< 
(29) 


We shall solve for a reliability function R(t) analogous to that given in Eq. (11) for the 
case of strictly periodic replacements. Let us assume that replacements can be made instan- 
taneously. In contrast to the situation treated in the last section, the number of replacements 
made in time t is a random variable. We shall show that R(t) is the solution of an integral 
equation of the kind encountered in renewal theory. 

The event "a system survives for a time t or greater" is the conjunction of events "the 
system survives for a time t or greater with j replacements having been made," where 
j= 1,2,3,... . Hence R(t) can be written, 





REPLACEMENT OF MACHINERY 


t 
a(t) = Rt) F(t) - | d F(z) R(7) R(t-7) F(t- 7) 
0 


t t 
| dF (7) { d F(T) -7,) R(7,) R(T9- 7) R(t - T) F(t - T9) 
0 7 


t t t 
- [ dF(7,) | d F(T) -7;) i dF (7, - 75) R(7,) R(T. - 7) R(T3- To) R(t- 73) F(t- 73) 
0 7 To 


We can simplify Eq. (30) by introducing two auxiliary functions defined by 


A (t)= R(t) F(t), 
(31) 


t 
w (t) = -F R(1) dF(1). 
0 


The function 0(t) is the probability that neither a failure nor a replacement occur in time t. 
The function W(t) is the probability that at least one replacement be made in time t. Equation 
(30) can be rewritten, by using Eq. (31), 


t 
R(t) = a+ [ Q (t- 7) dw (7) 
0 


t t-T 
+f dy (r,) { 1 ay (To - 74) 2(t- 79) 
0 0 


t t-T t-T 
+f aviry) [ avrg ey | 2 aw (11-79) 2 (t-73) 
0 0 


This expression for R(t) is seen to be the Neumann series, defining the solution of 


t 
(33) R()=act+ [ R(t- 7) dwW(7). 
0 


If R(t) is continuous and F(t) has.a derivative -g(t), then w(t) is differentiable. We shall let 
its derivative be denoted by @ (t). Under these conditions Eq. (33) is equivalent to 
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t 
R (t)= 2 (t) + R(t - r)$(r) dr. 
0 


Although it is seldom possible to find an explicit solution to either of Eqs. (33) and (34), 
we may still find the moments of R(t) by differentiating its Laplace transform at the origin. 
Let us define the following transforms: 


@ @ 
R*(s) = | eSt arity, at(s)= [ eSt gait), 
0 0 


@ 
v*(s)= [ eSt ay (t). 
0 


The Laplace transform of Eq. (33) yields for R*(s) the expression 


n* (s) 
(36) R* s) = —————__ + 
ts) 1 - y*(s) 

The moments of R(t) can now be recovered in terms of derivatives of R*(s): 


d™R* (s) 


(37) <t™> = (-1)" - 
ds s = 0+ 


In the simpler case where there exists the function ¢ (t) we may define another set of trans- 
forms: 


re) * 00 
Ry(8)= R(t) e"St at, ay(s)= [ a(t) eS" at, 
0 0 


i> 2] 
$*(s) = ¢(t) est at. 
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In this case we have 


(39) Rj (s) = 
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The moments may be written,*+ 


<t® = - { 


(+ 9) ie) 1 
t" aR(t)= n { t"-* R(t) dt 
0 


0 


* 
d®-1R. (s) 
= ( 1)" n A 
ds®-1 


* 
When the function R, (s) is used, computation of the first moment is particularly simple. As 
an example of the use of Eq. (40), let us consider the reliability function defined by Eq. (15), 
with F(t) equal to exp(- t). The functions 9, (s) and ¢ *(s) are given by 


1 


. - j+1 /n : 
a(s)= PCD G) sere 


*(s) = 8 a4 (s) ; 


For the first moment we find 


* 
ay (0) 
1-6 9; (0) 


<t>= 


n 
FL +48) 


Details of the integration are given in the Appendix. Several curves of ®<t>are plotted in 
Figure 4 as a function of a/f. 





4In deriving the second integralfrom the first, use was made of the fact that if <t™> is finite 
then lim R(t) t2? = 0. This is easily seen from the relations 
t-@ 
00 


@ 
t™ R(t) = - t? ake <- | x" dR (x). 
t t 


The existence of a finite <t™> insures that the last integral tends to zero for t ~@. 
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Figure 4 
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APPENDIX 


Derivation of Eq. (42) 


* 
Equation (41) indicates that om (0) is given by 


7, #-¢ aet mi 
(a1) 2 (= 2-1" (5) 


B+aj- 
We will now evaluate this sum. Consider a function J(u) given by 


, j-1 
T(u) = > (- 13+ 1 (?) eons 
j=l , 


=yP-t er kD 


It is easily seen that 


* rl 
Q, (0) = J(u) du 
1 

i 


1 
- { uf-1 ee du. 


1 
a 


We change the variable of integration by setting 


(A4) 


Equation (A3) then becomes 


11/1 @ 
(A5) a; (0) “a - Vv (1 - v)@ dv. 
0 


However, the last integral is just the Beta function [8], i.e., 
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* 1 1 
2 (0) =5 “+s p(4 








where the second line is obtained from the first by using the well-known expression for the 
Beta function in terms of Gamma functions. Hence, we have the expression 


which was to be proved. 
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NEWS AND MEMORANDA 


Readers are invited to submit to the Managing Editor items of general 
interest in the field of logistics. 


Dr. I. Heller, while on a leave of absence from the George Washington University 
logistics Research Project, has accepted a position as Mathematical Consultant with the Navy 
Management Office. 


The establishment of the Army Logistics Management Center at Fort Lee, Virginia, 
has been announced. The new Center represents an expansion of the Army Supply Management 
Course established in 1954 and will, in addition to the present course, offer four new manage- 
ment courses in the specific fields of Procurement, Requirements, Storage and Distribution, 
and Maintenance, with a basis orientation in the over-all supply system. 





The Secretary of the Navy has authorized the establishment of a new supply demand 
control point, the U. S. Navy Training Device Supply Office, effective 1 July 1956. The new 
SDCP is located at the U. S. Naval Training Device Center, Port Washington, New York, and 
will exercise inventory management over repair parts required to support all training devices, 
and over minor training devices, aids, and components as specifically assigned. The TDSO, 
commissioned on 2 July 1956, will supersede the Program Support Department of the U. S. 
Naval Training Device Center as the SDCP for cognizance symbol T material. The Center will 
continue to exercise technical control over this material. 





THE INDUSTRIAL COLLEGE OF THE ARMED FORCES 


The Army Industrial College was first established in 1924. In 1946 it was reconstituted 
as the Industrial College of the Armed Forces, under the supervision of the Joint Chiefs of Staff. 
The College, located at Fort McNair in Washington, D. C., conducts three closely related courses. 
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One of these is a 10-month resident course for senior military officers and civilian 
employees of the Federal Government. 


The second is a correspondence course called "Emergency Management of the Nationa] 
Economy." This course is open to qualified officers, both Regular and Reserve, and to civil- 
ians who have attended the National Resources Conference. The correspondence course 
normally takes a year to finish. Reserve officers receive 48 points toward retention, pro- 
motion, and retirement if they complete the course. 


The third course is the National Resources Conference, which seeks to foster a better 
understanding of the many interrelated problems associated with national security, a stronger 
appreciation of the inseparable nature of the civilian-military team, and a clear recognition of 
the responsibilities and capabilities of each half of the team. It seeks to encourage application 
of the highest level of skill, intelligence, and imagination to the complex problems which must 
be solved if our heritage of freedom is to be secure. 


The Industrial College is unique. It is the only military school which is principally 
concerned with the civilian economy and civilian-military relationships. All three Industrial 
College courses are based on the concept that, under our democratic form of government, over- 
all control of the military operations must be supplied by a civilian-controlled industrial com- 
plex; that moral, social, political, industrial psychological, and economic factors have a pro- 


found influence on military decisions. Both halves of the civilian-military team must understan 
these interrelated factors. 


The Industrial College is dedicated to the achievement of closer cooperation, more 


harmony, better understanding, greater sympathy, and smoother functioning of the civilian- 
military defense team. 


The College is neither a planning nor an operating agency. It is an educational institution 
whose purpose is to expand the nucleus of informed and thinking individuals who may be called 
upon to play important roles in the planning and managing of our national defense. Members of 
the College are not spokesmen for the Department of Defense or the military services. Their 
mission is to present some of the complex national security problems, many of which.are cur- 
rent and require solution now. They present no "canned" solutions, but rather some of the 
factors which must be considered in arriving at solutions. 
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THE SIGNAL CORPS: THE EMERGENCY. By Dulany Terrett, Office of the Chief of 
Military History, Department of the Army, Superintendent of Documents, Government Printing 
Office, Washington 25, D. C., 354 pp. 


"The more mobile an armed force becomes, the more rugged the terrain it encounters, 
or the more widely the force is deployed, the greater becomes the difficulty of securing and 
maintaining rapid, completely linked communications,"' General Orlando Ward points out in his 
Foreword to THE EMERGENCY, the first volume in the history of the Signal Corps in the Tech- 
nical Services subseries and thirty-second in the series UNITED STATES ARMY IN WORLD 
WAR II. He describes the Signal Corps as the agency which has, as part of its mission, the 
responsibility of developing, procuring, and furnishing signal equipment—no mean task in an 
age of swift and startling progress in electronics. 


The development of radar, second only to nuclear fission as the greatest scientific 
advance of the war, is the most important story in THE EMERGENCY. Here, for the first time, 
is given a complete account of the Signal Corps' part in the development of army radar. The 
volume also covers the development of frequency modulation and describes how, thanks to the 
Signal Corps, the U. S. armored force was the first to have FM, which all but revolutionized the 
use of tanks in World War II. The advance of crystal control is an episode linked with the FM 
story. 


The close relationship between the communications industry and the Signal Corps in the 
manufacturing of equipment and in the selection and training of officers and signalmen is 
described, as well as the wider but very much weaker relationship between the United States 
and its allies, especially the British. In connection with the latter, the mission of Sir Henry 
Tizard and other electronic scientists and physicists to the United States in 1940 is covered. 


Part One of THE EMERGENCY narrates the history of the Signal Corps from its earliest 
days up to World War II. It describes the Corps in World War I, postwar curtailment and peace- 
time procurement planning, the scope of the Signal Corps during this period, and early research 
in wire, radio, and radar. Part Two continues the story up to the eve of Pearl Harbor and 
covers the rapid expansion of the Signal Corps during the Limited and Unlimited Emergencies. 
The remarkable developments in the electronics field are described, as are the various Serv- 
ices—Pigeon, Photographic, Tactical and Administrative, and Supply—of the Corps. 
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A thirty-six-page appendix of World War II Signal Corps equipment will prove of 
interest to the communications expert and amateur alike. 


THE TRANSPORTATION CORPS: MOVEMENTS, TRAINING, AND SUPPLY. By 
Chester Wardlow, Office of the Chief of Military History, Department of the Army, Super- 
intendent of Documents, Government Printing Office, Washington 25, D. C. 


In this volume of the series U. S. ARMY IN WORLD WAR II, emphasis is given to the 
fact that, in a great conflict fought on foreign soil, victory depends on how many soldiers and 
how much materiel can be moved to the oversea commands and on the timeliness of their 
arrival. 


It becomes apparent from this history that supplying the oversea theaters involved far 
more than the transportation of huge tonnages. (The Army shipped about 132,000,000 measure- 
ment tons from U. S. ports during the war.) The task was an intricate one of coordination and 
synchronization. The right proportion had to be maintained between troopship and cargo ship 
capacities. Care had to be taken to ship to the theaters the types of troops and the items of 
supply that they needed and to deliver them when they were wanted and at ports where they 
could be accommodated. To assemble and ship related supply items received from many dif- 
ferent depots and factories required careful identification and sorting at the ports of embar- 
kation. The great variety in weight and size of military supplies made it impossible to get 
well-balanced cargoes and efficient use of cargo space, and the author observes that, ''. . . the 
tendency of the commanders of active theaters, fighting thousands of miles from their main 
source of supply, was to place heavy requisitions in order to safeguard themselves against 
unforeseen supply requirements or an interruption in the line of communications . . . some 
oversea commanders accumulated larger stocks than they could properly administer . . 

This situation, coupled with the willingness of the theaters to detain vessels for long periods 
and use them as floating warehouses, led the Chief of Transportation to assert that the accu- 
mulation of unnecessarily large stocks overseas was one of the logistical mistakes of the war." 


The movement of troops and equipment within the United States is an integral part of 
the transportation story told in this volume, and the author describes how the Military Trans- 
portation Section of the Association of American Railroads, set up early in the emergency, 
functioned side by side with the Chief of Transportation's staff in working out efficient routings, 
conserving equipment, and maintaining schedules. 


The relationship of the Transportation Corps with the Office of Defense Transportation, 
the wartime agency charged with the control of inland transportation, and with the War Shipping 
Administration, which controlled the bulk of merchant shipping, comes in for careful scrutiny. 
The book also describes the Corps' struggle to resist attempts by powerful elements in the 
Army Service Forces to take over functions that had been assigned to the Transportation 
Corps. 


While the movement of troops and materiel was the distinctive function of the Chief of 
Transportation, he had other exacting responsibilities. Such tasks as training troops to operate 
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ports, railroads, small boats, and amphibious trucks and to handle other transportation activi- 
ties in the overseas commands became increasingly heavy as the war progressed. Designing 
and procuring the many types of small vessels required by the Army and procuring locomotives 
and rolling stock for the military railroads and for lend-lease were additional major concerns 
of the Chief of Transportation. A full account of these aspects of the transportation job appears 
in the volume. 


MEASURING BUSINESS CHANGES. By Richard M. Snyder, John Wiley and Sons, Inc., 
New York, 1955. 


The purpose of this book is to serve as a handbook of significant business indicators or 
indexes, it consists of a general description of the indicators, plus the details included in them. 
The indicators discussed are: (1) national income and product, (2) population, (3) labor: 
employment and earnings, (4) commodity prices, (5) production and business activity, (6) con- 
struction activity and costs, (7) trade, (8) financial activity, and (9) stock prices. The indica- 
tors described are primarily those published by government agencies. 


MAINSPRINGS OF THE GERMAN REVIVAL. By Henry C. Wallich, Yale University 
Press, New Haven, 1955. 


Professor Wallich views the factors which have contributed to the economic recovery 
and growth of West Germany since World War II. He considers the domestic economic and 
political policies, as well as the effects of the economic and political events which occurred in 
the outside world during the post-war period. One of the most interesting aspects of the work 
is the author's account of the economic recovery and growth through the existence of a rela- 
tively free market policy and the government's "tight money" and balanced budget programs, 
while most other nations were proceeding in other directions to foster economic development. 


THE ECONOMICS OF NATIONAL AND MILITARY STRATEGY. By Robert E. Kuenne, 
United States Naval War College, 1956. 


This book is a compilation of seventeen lectures delivered by the author in 1956 in imple- 
menting the economic portion of National Security Studies at the Naval War College. The 


lectures are,therefore, structured, with the intent of developing a framework for the analysis of 
economic problems arising from the planning and implementation of foreign policy. As a 


result, the book has the aims of a textbook in a field in which few attempts have been made to 
bring the tools of modern economic analysis to bear. 


In Part I-- termed the “Economic Factors” --a framework is developed for the appraisal 
of an economy’s ability to support an ambitious foreign strategy. Ten such “variables” in the 


strategic equation are isolated, an attempt is made to indicate their role and relative strength in 
the equation, and five of them are set apart for more detailed study within the context of the 


Present world power struggle. These latter factors are the level of national product, the rate 
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of growth of national product, the allocation of resources, the spatial and temporal structures y 

an economy, and the nature and degree of external economic linkages of national economies, 
Part II constitutes an attempt by the author to judge the present levels of these five 

variables for the United States economy, and Part III is an attempt to answer the same questions 


for the Soviet economy, Having presented in Part I a framework relating strategic variables of 
an economic character to the development of strategy, Parts II and III try to “plug in” the “valyg;’ 


of these five most important variables into the strategic equation in order to obtain some idea of 
present and future capabilities of the two major power blocs. 

The author attempts to demonstrate the manner in which these five variables can be 
manipulated in periods of cold, limited, and total war in order to achieve desired strategic 


capabilities. This analysis is contained in Part IV--“Economic and Strategic Thinking” -- and 
includes the discussions of tax policy, deficit speding, mobilization problems, the use of input- 


output and linear programming techniques in such manipulation, economic warfare in the con- 
ventional sense, and foreign aid as a strategic instrument. 


Lastly, in Part V, problems of economic development and their meanings in terms of 
shaping strategy are discussed in specific terms of the conflicting economic programs of India 
and China. _ 

The author is assistant professor economics, Princeton University, who served as con- 
sultant to the Naval War College in 1955 and 1956. 


MODERN MATHEMATICS FOR THE ENGINEER; Edited by Edwin F. Beckenbach, 
McGraw-Hill Book Company, Inc., 1957. 


This book gives a good quick look at what is new and exciting in the world of highbrow 
mathematics doing lowbrow chores. Opening with nonlinear mechanics and closing with high- 


speed computers, it touches all the high points along the way — such as calculus of variations, 
partial differential equations, theories of prediction, games, operations research, and computing 


methods. There are nineteen papers arranged in three groups: Mathematical Models, Probabil- 
istic Problems, and Computational Considerations. Each is written by an authority (frequently 


the authority) on the subject covered. This book will be valuable to the engineer who is opening 
the literature on a new subject. The papers themselves plus the excellent reference lists will 


take him at once into what is new and important. With the compression required to cover this 
wide range in about 500 pages, it is clear that the treatment is not exhaustive, but even so, it is 


likely that many readers will find new material and ideas in subjects they have thought they 
knew well. 


FASTER, FASTER. By W. J. Eckert and Rebecca Jones, McGraw-Hill Book Co., Inc., 
1957. 


“A simple description of a Giant Electronic Calculator and the problems it solves.” 


Definitely recommended reading for everyone who has an interest in computing machines. 
The authors do not tell you how to design and build a computer but they tell you why computers 
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are designed and built as they are. They do not give a course in programming and coding but 
they describe how a problem is programmed and coded. And most important, they devote one 
chapter to the type of problems appropriate to computer calculation. For too long the subject 
of computers has been covered by monographs either too technical or too elementary for the 
average reader. In Faster, Faster the authors have hit a nice balance in their presentation of 
material and this book should be on the bookshelf of every mathematician. 
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